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NECROLOGY 


W. A. GODDIJN 
(April 9, 1884-February 22, 1960) 


To few of us it is allowed to play an important part in a department 
of science outside our official field and outside the field of our 
university teaching. GODDIJN was one of them. A pharmaceutical 
chemist by profession, his true interest went out to the application 
of genetical methods in taxonomy. 

Wouter GoppiJN studied pharmacy at Leiden. After his gradu- 
ation in 1909 he became an assistant at the Government Herbarium 
at Leiden, and in 1914 accepted the appointment of Curator. In 
1926 he obtained the degree of Ph. D. with a thesis concerning a 
pharmaceutical plant (Hyoscyamus niger). In 1934 he was appointed 
Associate Professor and in 1939 Full Professor to teach pharma- 
cography, galenic pharmacy and the methods of dispensing. 

Although he performed all his duties with that painstaking de- 
votion which was such a characteristic feature of his personality, 
GOoppIJN'’s heart was not in pharmacy, but in botanical taxonomy 
and genetics. 

He made numerous genetical experiments within the genera Vrola, 
Gewm, and particularly Scrophularia. In Scrophularia he obtained 
many years ago probably the first example of a truebreeding artificial 
allotetraploid. He never published this result, because the director of 
the Government Herbarium at that time did not consider his research 
completed. In the company of J. P. Lorsy GoppijN travelled in 
South Africa to study hybrids in the native flora and hybrids in the 
human population (Genetica, 10, 1928, 315 pp. + 11 pl). He also 
studied genetical problems in the floras of Switzerland and Italy. 
His “spare time” was devoted to a self-imposed and thankless task, 
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invisible for outsiders: from 1930 to 1959 he was the general editorial 
secretary of the Dutch genetical periodicals. 

GODDIJN did not like to be in the foreground; he never asked for 
appreciation. Nevertheless he earned the greatest esteem, especially 
of those who were in a position to closely follow his activities. His 
poor health did not impede his work; his apparently weak body 
housed a strong mind, forceful and pure, honest and upright. 

WOUTER GODDIJN was a brilliant scientist, an extremely modest 
personality and an uncommonly good character. The old saying: 
“All for love, nothing for reward”’ applies to him in its fullest sense. 


M. J. Sirks 
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THE SEPARATION OF EPISTATIC FROM ADDITIVE 
AND DOMINANCE VARIATION IN GENERATION MEANS. II 


by 
B. 1. HAYMAN 
Applied Mathematics Laboratory, D.S.I.R., Private Bag, 
Christchurch, New Zealand 


(Received for publication January 23, 1960) 


INTRODUCTION 


The previous article (HAyYMAN, 1958) pointed out that additive 
and dominance variation cannot be uniquely specified when the 
genes concerned interact with each other. Consequently, genetic 
variation cannot be partitioned uniquely into additive, dominance 
and epistatic components and the relative contributions of these 
modes of gene action to various genetic phenomena such as heterosis 
cannot be ascertained. This article supplies an alternative approach 
to partitioning problems which is a development of a suggestion in 
the previous article and in an article on diallel crosses (HAYMAN, 
1957). It also demonstrates how misleading the partitioning method 
can be. 


PARTITIONING THEORY 


In the previous article the difficulty in specifying additive and 
dominance variation in the presence of epistasis was brought out 
in two ways. Since this difficulty is fundamental to our discussion 
we present further demonstrations of it here. In the case of two 
genes controlling a quantitative character the model of HAYMAN 
and MATHER (1955) gives the expectations of the nine possible 
phenotypes in terms of nine parameters, m the mean, da and dy 
measuring additive effects of genes a and b, ha and /, measuring 


— 133 — 


134 B.I. HAYMAN, THE SEPARATION OF EPISTATIC FROM ADDITIVE 


dominance effects and #av, fav, Îoa and lap measuring interactions 
of da and he with dy and Aj. In this model the additive effect of 
gene a, defined as half the difference in effect of its two homozygotes, 
can take three values, 

da + dav — HJab 

da + Fav 

da — lab — Hiab 
depending on whether the other locus is occupied by BB, Bb or bb. 
The average of these three quantities in the proportions 1 : 2: 1 
is da; in other words, da in this model of gene action is the additive 
effect of gene a against an Fz background population at locus &. 
An F2 background is, of course, only one of many that may be 
arbitrarily chosen, some being more appropriate in certain genetic 
experiments than others, but with all other background populations 
the additive effect of gene a differs from the d, defined above by 
multiples of its interactions, Zap and jan, with gene b. 

In the case of several genes controlling the character in question 
HAYMAN (1954) has shown that if mm, da, ha, lap, Jav and lap are 
defined against a background population containing AA, Aa and aa 
in the proportions pa : qa : Ya and if m, da, ha, Zap, dak and Zap are 
defined against background proportions pa : da : Fa then 
mM = MJ pe («ada je Baha) sie Ai (aaartab st «ab vab el «pPaîva de 

À gen + Bafvlav) 


da — da + 2 (artan + Êvîav) == Omlôaa 
b 

ha = ha + X (avfva + Blan) — Omlopa 
b 


Tab — tab == O2mfdaadar 
Jar = jev — 02m lop 
lav — lab = O2ml0B0Bs 
where «a = Pa — Fa — Pa + 7a 

Pa — Ja — da 


Interactions between three or more genes have been omitted but 
it is clear how the formulae may be generalised to include these. 
The three earlier values for the additive parameter are special cases 
of this general formula for dà when (ba, qa, ra) has the value (4, 4) 
and (Pa, Ja, Fa) takes the successive values (1, 0,0). (O, 1,0) and 
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(0, 9, 1). J. H. VAN DER VEEN (1959) has given some special cases of 
these formulae. 

We will restrict our discussion of experiments to the cross, F3, 
between two inbred lines, P and P’, and descendant families obtained 
by selfing and backcrossing, etc. The only statistics considered in 
this article will be family and generation means. 

The components of means are obtained as follows. Let the para- 
meter Ô4 take the value of positive or negative unity according as 
parent P contains the positive or negative homozygote of gene a. 
Let (bag, Jag, ag) be the proportions of AA, Aa and aa in generation 
g descended from the two inbred lines. Table 1 contains the values 
of Pag — Yag And qag in various generations, the former being a 


TABLE 1. Gene frequencies in genevations 
descended from two inbved lines 


Generation Pag — Yag Jag 
| 

P Oa 0 
P’ —0a 0 
Fi 0 Î 
Fa 0 5 
F3 0 4 
Fa 0 E3 
B 30. È 
B! — 3 
BS Za or 
BS’ — 30. ï 


multiple of 9, and the latter numerical. If (ba, ga, Ya) is the common 
background population against which the genetic parameters are 
defined in all generations then in the absence of linkage the mean 


of a generation is 
Mg S= M Se DA (aagda in Bagha) de 
a p n 2 
+ N (aagxvgtav zn avagPvglab ae apgPaglva zijn BagBvalav) 


a<b 


where aag — Pag — Pag — Pa + Ya 
Bag — dag — Ya 


If the background population is one of the generations under con- 
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sideration, or any average of these generations such as the midparent, 
we may write 


Ag S= caglÔa 
Êg — Pag for all loci 
and d => Oda 
h == DE ha 
1 = D) OaOvtav 
a<b 

j == 4 X (Oafav + Ôofva) 
Ll = DD lab 


so that mg — m + agd + Byh + agú + Zagpai + Bot. 

This also establishes HAYMAN's (1958) unproved formula for the 
generation means of this experiment in terms of the components 
m, d, h,i, j and /. It can obviously be extended to-include higher 
order epistasis. 

HAYMAN’s (1958) rule for writing down the expectations of gener- 
ation means in terms of different background populations follows 
from the fact that «4 and 9, are zero in the expectation of the mean 
of the background population. This expectation is taken to be 7 
alone and all other expectations have d and A terms added with 
appropriate coefficients, either from the formulae for ag and fy or 
from the well-known relationships between the generations. Epistatic 
terms to the order desired are added according to the above formula 
for mg or its generalisation. 

Three points emerge from this discussion. The first is that no 
background population has an advantage over the others either on 
genetical grounds or on grounds of algebraic simplicity. The second 
point is that change of parameter from one background to another 
is determined by interactions involving that parameter. The value 
of ha is only affected by jpa and lap and not by tav or jav. Conse- 
quently, if interactions involving more than, say, k genes are absent 
then interactions involving k genes are determined uniquely and 
independently of the background population. Thus d and % are 
defined uniquely in the absence of any epistasis and so are ú, j and 
lin the absence of interactions between three or more genes. The 
third point is that these definitional differences are of the same 
order of magnitude as the interactions involved. In other words, a 
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parameter is indeterminate to the extent that its interactions exist 
and, in particular, if epistasis is as important as additive and domi- 
nance variation then the latter two quantities have little meaning 
as distinct entities. 

Another and shorter illustration of the difficulty of defining 
additive and dominance variation in the presence of epistasis is 
the analogy between epistatic and non-epistatic variation on the 
one hand and non-epistatic variation and the mean on the other. 
In a non-epistatic system the expectations of two inbred lines and 
their hybrid may be written 


P =mtd 
PS md 
Fi =mh 

or P =md-d—h 
P' =m—d—h 
Fi =m 


depending on whether the mean 1 is situated at the mid-parent 
(background population) or at the Fi (background population). 
Both these models give unique estimates of d and Jh as 

d=HP—P) 

h=Fi-HP+AP) 
but their estimates of 1» are different. Indeed the mean may be 
situated anywhere without altering the estimators of d and Jh. Evi- 
dently the mean m has no intrinsic genetic meaning and is merely 
a statistical complement to the genetically meaningful quantities d 
and A. If, in this last statement, we replace m by m, d and h, and d and J 
by #, j and / we obtain a statement valid when epistasis is present. 

To sum up, to attempt to partition the variation between the 

generation means amongst components m, d, h, i, j and ! is a meaning- 
less exercise if the epistatic components are important. This con- 
clusion may seem obvious but, before discussing an alternative 
approach to the partitioning problem, we will reinforce our con- 
clusion with a numerical example of the misleading consequences of 
partitioning against an arbitrary background population. 


EXAMPLE 


The problem is the partitioning of heterosis and the example is 
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the cross between lines 1 and 4 of Nicotiana rustica grown in a diallel 
cross of eight lines in Birmingham, England, in 1952. Various papers 
by JINKs and HAYMAN have discussed these crosses but the data of 
interest here are in a recent article by JINKs and JoNEs (1958, 
table 6). 

Heterosis is the situation of the Fj outside the range of its parents 
and may be measured by 


F1 — max (P, P’) 
if positive, as in our example, or by 
F1 — min (P, P') 


if negative. Heterosis may be partitioned into components against 
any background population. We choose four backgrounds, the 
midparent, the parent P, the cross Fi and its selfed progeny Fo, 
the latter being used implicitly by JiNKs and JoNes. If m, d, h, #, 
j and / are the mean, the additive, the dominance and the epistatic 
components of family or generation means defined against the first 
background population and if one, two and three primes indicate 
components defined successively against the other background popu- 
lations then 


m=mdtdhi m=mthtl m'=mtihtH 
LD d' =d + 2 d" =d +4 
h =h +2 hhh GA KEE, 


The experiment was carried to F9 and backeross generations which 
provide an exact fit to these components and no test for higher 
order epistases, which are therefore ignored. The epistatic com- 
ponents, £, j and /, are determined uniquely but not the mean nor 
the additive nor dominance components. Heterosis may be partition- 
ed as 
Fi -P=d thi 1 

=d HW Hil 

=d AW il 

NN sie nh zedig 
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The values of JINKs and Jones for the cross 1 X 4 were 


Eb 
d' = _ 2.05 
h'—= _ 4,25 

i= —16.20 
1 == —6.85 
= 15.32 


From these have been derived the alternative partitions of heterosis 
in table 2. The significance levels of the components are determined 
by comparison with standard errors averaged over all crosses of the 
diallel and so may vary a little from the significance levels of JINKs 


TABLE 2. Components of heterosis in cross T X 4 against different 
background populations 


Background population 
Component 
Sym. hom. 2 Fi F3 
Additive —8.90t** 23.50 4.80 —2.05 
Dominance — 11.07 — 24.77 191574 4.25 
Total non-epistatic —19.97 —1.27 24.37** 2.20 
add. x add. 16.20* —16.20* 16.20* 16.20* 
Epistatic add. x dom. — 13:70 —13.70* —6.85* 
hi dom. X dom. |__ 15.32 1532 | —15.32 n 
total Sl o2 12.82 — 12.82 995 
Heterotic U heee iid IES (OSE 


and Jones. However, the situation is clear from the means alone. 
The last column of figures provides JINKS and JONES’ picture of the 
genetic situation in this cross — heterosis apparently arises from % 
interactions assisted by a small dominance contribution and opposed 
by small j interactions. The other columns of figures present different 
pictures. In the first column epistasis is still the major factor in the 
heterosis but it comprises two equally strong components, # and /, 
and is opposed by both the additive and the dominance components. 
In the second column additive and dominance effects just about 
cancel out leaving a balance of all three epistatic components as 
the source of heterosis. In striking contrast, the third partition 
contains the same epistatic components as the second, but reversed 
in sign so as to oppose the heterosis which is here largely due to 
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dominance. Bringing the significance levels into the picture does 
not alter the conclusion that partitions such as those in table 2 can 
provide no information about the relative importance of additivity, 
dominance and epistasis to heterosis. 

Another approach to this problem of partitioning heterosis which 
avoids the consideration of varying background populations is to 
form the regression of heterosis on its possible and likely components. 
Since the additive and dominance components are not well-defined 
in the presence of epistasis we calculate the regression of heterosis 
on the three epistatic components alone. This cannot supply a 
complete predictive equation for heterosis but it should indicate, 
by evaluation over a number of experiments, whether there is a 
consistent relation between heterosis and the well-defined epistatic 
components. Such a relation would be evidence for a general par- 
tition of heterosis into additive, dominance and epistatic components. 

For this test all 28 crosses of JINKS' diallel experiment were used 
and regressions of F, — max (P, P’) on #, j and / calculated for the 
two characters, plant height and flowering date, recorded in repe- 
titions of the experiment in 1952 and 1953. Epistasis was present in 
all four cases (P < .001 of no epistasis). Heterosis only occurred in 
height, and then not in all crosses, but in the majority of crosses in 
both characters the F, deviated significantly from the greater parent. 
Continuity of gene action indicates that any model capable of 
explaining heterosis should explain all these deviations whatever 
their sign so that we have here four estimates of such a model. The 
results are expressed in terms of regression equations and of the 
multiple correlation coefficients between observed and expected 
heterosis which test the significance of the regressions. 


Height, 1952: r == 0.594* 
heterosis — 1.603 — 0.365* — 0.152j — 0.319** 
Height, 1953: LOA 


heterosis — 1.892 + 0.058i — 1.426***j — 0.072/ 
Flowering date, 1952: 7 = 0.173 
deviation —= —4.412 — 0.116i — 0.123j — 0.0671 
Flowering date, 1953: 7 —= 0.521 
deviation — —14.547 + 0.2361 + LO16**j + 0.221 
The asterisks denote significance at 5%, 1% and 0.1% probability. 
The final correlation is nearly significant at 5%. 
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There is no consistency in these results. In height, # and / inter- 
actions determine heterosis in one year and j interactions in the 
other year. In flowering date one regression gives no evidence of 
any relation between deviation of F1 from parent and epistasis 
although both the deviation and the epistasis are highly significant 
in many crosses. The other regression of flowering date suggests 
that j interactions may be important although the overall regression 
is barely significant. It is clear that this approach supports our 
previous conclusion that heterosis cannot be partitioned uniquely 
into additive, dominance and epistatic compounds such as we have 
been using. 


LEAST SQUARES THEORY 


In the partitioning of family means described above only com- 
ponents of highest order had a useful interpretation; lower order 
components were indeterminate because of their interactions. We 
will show how this indeterminancy may be measured and how the 
components may be redefined to restore their usefulness when their 
interactions are not large. 

We can classify the simpler genetic systems into three groups. 


Group |. No significant epistasis. Significant additivity and/or domi- 
nance. 

Group 2. Significant epistasis. Additivity and/or dominance signifi- 
cantly more important than epistasis. 

Group 3. Significant epistasis, additivity and/or dominance, with 
epistasis at least as important as the other modes of gene 
action. 


Group 1 presents no problem of estimation and is trivial as far 
as this article is concerned. Group 3 contains the genetical systems 
to which the earlier discussion of partitioning is appropriate. Little 
meaning can be attached to the concepts of additivity and dominance 
in such systems but the influence of epistasis will be discussed later. 

Group 2 is intermediate and merges with the other two groups 
at its extremes. Genetic systems in this group are characterised by 
significant epistasis which is nevertheless negligible in importance, 
compared with additivity and dominance, in explaining variation be- 
tween the generation means. They are akin to the agricultural trials 
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which by their statistical sensitivity detect significant differences 
between treatments, or between varieties, which are too small to 
be of economic importance. In these genetic systems it is best for 
a start to ignore the epistasis so as to obtain unique estimates of 
additivity and dominance alone and to incorporate the epistatic 
indeterminancy, together with the usual errors of environment, 
measurement and genetical sampling, in their standard deviations. 

We follow the approach of HAYMAN (1957, 1958) and use the same 
estimators as in group |. We initially ignore the significant but 
relatively trivial epistasis and estimate d and h by least squares 
from the equations. 


Mg = M + agd + Boh. 


The only difficulty is to determine the standard deviations of the 
estimators. In group l these are calculated from the error mean 
square in the unweighted least squares analysis or from the weights 
in the weighted least squares analysis. In group 2 the unweighted 
analysis exhibits a hierarchy of mean squares for additivity and 
dominance, for epistasis and for error, the first being significantly 
greater than the second and the second than the third. Then the mean 
square for epistasis replaces the error mean square in determining 
the errors of the estimators. In the more usual weighted analysis 
the hierarchy consists of chi-squares for additivity and dominance 
and for epistasis, both being significant but the former being sig- 
nificantly greater than the latter. The standard deviations of the 
estimators of additivity and dominance are functions of the weights 
and it is suggested, at least tentatively, and by analogy with the 
unweighted analysis, that when epistasis is present these functions 
be multiplied by xf-* where f is the number of degrees of freedom 
of the 7? indicating significant epistasis. 

Statements about these unique additive and dominance components 
should be reasonably correct as long as their significance is adjudged 
on the corrected standard deviations. When epistasis is absent the 
genetic system moves into group 1 and the statements are certainly 
correct. When epistasis is as important as additivity and dominance 
the genetic system moves into group 3 and the standard deviations 
correctly ensure that no significant statement can be made about 
additivity and dominance. 
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Formulae have been published for estimating additive and domi- 
nance components of means in two types of genetic experiment. 
One (HAYMAN, 1957) concerns parental, Fi and F3 families of a 
diallel cross of inbred lines and the other (HAYMAN, 1958) concerns 
parental, F3, Fa and back€ross families of a single cross between 
inbred lines. The accompanying standard deviations should be 
corrected upward as suggested here. 

The Nicotiana rustica cross cited earlier is a genetic system in 
group 2. Table 3 contains the chi-square analysis of the variation 
between family means. Epistatic variation is significant and additive 
and dominance variation is significantly greater. When components 
m,d and Jh are fitted 
L 2.71 
LL 4,46 


and these two components account for the major part of the genetic 
variation between means. The variance formulae have been multiplied 
by 4.88, the mean square for epistasis in table 3. 


TABLE 3. Analysis of chi-square for cvoss IT X 4 


Variation | eeens | Degrees of | Mean | fe 
category | freedom square 
Add. + Dom. | 136:35*2* | 2 | 68.18 | 1397 
Epistatic | 14.63** 3 4.88 
Total genetic | 15098222 | 5 | | 


The answer to our problem of the source of heterosis is now 
clearer. On the evidence of parental, F, Fz and backcross means, 
heterosis in this cross can be ascribed in the main to dominance. 
Whether this is overdominance, or repulsion association in the 
parents, cannot be ascertained without reference to segregational 
variation within families. This is beyond the scope of this article. 

The next step in the analysis of genetic systems in group 2 is 
to estimate m, d and Jh defined against some background population, 
together with #, 7 and /. The differences between observation and 
expectation in the previous model ignoring epistasis are epistatic 
in nature and their expectations in the epistatic model may be 
expressed in terms of these i, j and /. These expressions reveal the 
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direction of influence of epistasis on the means and the relative 
importance of 4, j and /. 

For our illustrative Nicotiana rustica cross the general expectations 
of the deviations of the means have already been given by HAYMAN 
(1958, table 1B). For the F1 mean, for example, 


(FiFi 02200 Abil 


With values of —16.20 and 15.32 for # and / respectively this devi- 
ation is seen to be the resultant of opposing epistases. The numerical 


value is : 
Fr Ar lS3 5077 


giving some evidence that epistasis diminishes the F} mean. 

We can complete the solution of the heterosis problem in this 
cross. The difference between the observed heterosis and that ex- 
pected on ignoring epistasis is 


Frabosa 


and its expectation is 0.2592j. The expected heterosis is signifi- 
cantly greater than that observed and the difference is due entirely 
to 7 interactions. Epistasis, as far as it is active, decreases heterosis 
in this cross. This is contrary to the result of JiNKs and JONES 
(1958) given in the last column of our table 2. It is more like the 
partition against an F, background population in the adjacent 
column but there is no reason to prefer any partition either in 
table 2 or from amongst the other possible partitions. 

The analysis of genetic systems in group 3 formally follows the 
lines of that for group 2 but no genetic interpretation is placed 
upon the estimates of additivity and dominance. Estimates obtained 
by ignoring epistasis are discounted by their standard deviations 
and those obtained in conjunction with epistatic components are dis- 
counted by their indeterminate backgrounds. The investigation is 
confined to comparisons between #, j and / and to their contribution to 
the differences between observed means and expectations ignoring 
epistasis. 

In summary, the weighted analysis of the three groups of genetic 
systems proceeds as follows. In group 1, additive and dominance 
components are estimated together with their standard deviations by 
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weighted least squares. The analysis in group 2 is initially the same 
as in group | except that the standard deviations of the two genetic 
components must be corrected upwards. In addition, epistatic com- 
ponents are estimated and the deviations of observation from expec- 
tation on the model ignoring epistasis considered in terms of the 
epistatic components. In group 3 the comparison of observed and 
expected means seems to be the only source of information. This 
reveals the relative importance of the three types of epistasis but 
cannot relate them to additivity and dominance which, indeed, are 
not realisable concepts when epistasis is a major source of genetic 
variation. 

This analysis may be extended to genetic systems exhibiting 
epistasis between three or more genes. The basis is again the hier- 
archy of mean squares or chi-squares. This reveals whether any 
particular level of gene action is relatively independent of its inter- 
actions, and components should be estimated for the lowest level 
with this property. 


CORRECTION 


In the previous article (HAyYMAN, 1958) it was suggested that, 
when parameters were defined against an Fz background pouplation 
and estimated from P, P’, F,, F2, B and B’ families, the estimators 
of m, d, hi and j were unbiassed in the presence of linkage. VAN 
DER VEEN (1959) has pointed out that this is incorrect. In fact, 
only estimators derived from P, P’ and Fi (the non-segregating 
generations) are unbiassed by linkage. Examples are 


and any linear combinations of these quantities. 


SUMMARY 


The difficulty of specifying additive and dominance variation in 
the presence of epistasis is described. This difficulty can be over- 
come when the epistasis is at a relatively low level of intensity but 
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when epistasis is the major source of genetic variation no measure 
of additive or dominance variation is possible. 

Heterosis in a cross between two inbred lines of Nicotiana rustica 
was found to be due to dominance opposed to a small extent by 
epistasis. The misleading nature of another method of partitioning 
heterosis is pointed out. 
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THE APPLICATION, TO TETRAD-ANALYSIS-DATA FROM 
SACCHAROMYCES, OF PRINCIPLES FOR ESTABLISHING 
THE LINEAR ORDER OF GENETIC FACTORS *) 


by 


ERNEST E. SHULT and SHARON DESBOROUGH 
Southern Illinois University Carbondale, Illinois 


(Received for publication February 4, 1960) 
INTRODUCTION 


1. Remarks on Gene-Order 

The first criterion utilized to establish gene order was developed by 
STURTEVANT (1913, 1915). This criterion defined the correct gene order 
of three genes ABC, as (a) that arrangement which yields the lowest 
frequency of double recombination. For example, consider the three- 
point cross ABC /abc Q x abc/abc 3, with the following distribution of 


TaBre 1. Genotypic Frequency 


Class Frequency 

@% ABCHRTSON 
abc 

or; Abc 9% 
aBC 

C3 ABc 40% 
abC 

Ca AbC 25% 
aBe 


*) This work has been supported by research grants from the U.S. Public 
Health Service, C-4682 and the American Cancer Society, E 125A. 
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genotypic classes. By criterion (a), the correct gene order is B-A-—C, 
and class Cz is the double recombinant class. 

The coincidence index is commonly defined as Ca/(C3 + Ca) (Ca + Co). 
This criterion (which will hereafter be designated the Sturtevant Cri- 
terion) may be defined in two other ways which are logically equiva- 
lent to (a): The correct gene order is that permutation of A-B-C, 
which (b) minimizes the coincidence index, or (c) which minimizes the 
value of 

f(A-B) + f(B-C) — f(4-0), (1.1) 


where f(4-B) is the frequency of recombination between A and B. 
(Note that (1.1) is twice the observed number of progeny classified 
as double exchanges when the order is ABC). The equivalence with 
(c) follows directly from the formula of Trow (1916), for the correct 
order minimizes (1-2f(A-C)) — (1-2f(4-B))(1-2f(B-C)). 

The Sturtevant Criterion is one of the most basic concepts of modern 
genetics, and is fundamental to the construction of chromosome maps. 
Some writers have considered that this criterion proves that genes are 
in a linear order. This is impossible since the criterion defines order — 
that is, the “order”, B-A-C is the name of those conditions whereby 
the arrangement B-A-C leads to the recognition of (a), the lowest 
frequency of double recombination, (b), the lowest coincidence index, 
or (c), the lowest value of (1.1). Second, application of this criterion 
automatically imposes an order on axy three transmissable factors 
(unless, Ca = C3 < C4), whether they are physically ordered or not. 

The fact which appears to establish the linear order of genes as a 
physical reality is the consistency of those orders obtained from a large 
number of three-point tests involving overlapping combinations of 
genes. The degree of certainty in each case is a function of the number 
of independent checks of gene order. Thus, with three genes there is 
only one possible 3-point test. The total information of this test is 
exhausted in defining gene order so that there is no means of confir- 
mation. With four genes, one three-point test may establish the order 
A-B-C, a second test, the order B-C-D, while the remaining 3-point 
tests must then confirm the predicted orders A-C-D and 4-B-D. 
Thus, this system has two checks. There are 


0 checks for 3 genes, 
2 checks for 4 genes, 
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7 checks for 5 genes, 
16 checks for 6 genes, 
30 checks for 7 genes, 
50 checks for 8 genes, and in general, 
n 


3l(n — 3)! 


(n — 2) = Car (1 — 2) checks for n genes. (1.2) 


As can be seen, the certainty of the gene order rises quite rapidly as 
the number of genes increases. 

Although this procedure has worked successfully in most organisms 
(e.g. Drosophila, Maize, Neurospora, Aspergillus, mice, etc.), it has 
achieved renewed interest in its application to phage genetics. It 
should be clear from the above discussion that a single three-point 
analysis (as in table I) is not sufficient evidence for the existence of a 
physically linear gene order. 


2. Genetic Distance and Gene Order 

Another assumption, fundamental to the development of modern 
genetics is that recombination is due to an event, Y, occurring at a 
particular point on the linear continuum between the genes. If only 
two parental alternatives are available (as is usual), two events of 
this kind leave distal points in parental combination. Note that this 
notion already presumes the principle of order, discussed in section 1, 
for the hypothesis states that two genes will recombine only if the 
event Y, and the genes A and B are in the order A-Y-B. “Genetic 
Distance” is then defined as the average number of such events, 
occurring between A and B. It is not necessary that genetic distance 
be a linear function of physical distance since the density of events 
along the physical continuum need not be uniform. However, genetic 
distance must correspond with physical distance in one fundamental 
property: both must be additive measures. It is only by virtue of this 
property that such an index earns the name “distance”. This con- 
dition is fulfilled for “genetic distance’, as defined above, since the 
mean of a sum of random variables is the sum of their means, irre- 
spective of whether random variables are dependent or independent. 
That is, E(X Aat … 4 Ka) = E(X2) 4... + E(X) where 
each variable, X;, assumes values, O, 1, 2, 3, ..…, equal to the number 
of events, Y, occurring in its corresponding subregion, Ri, and E(X;) 
is the expected value (mean) of the random variable, X4. 
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In virus genetics, a DNA or RNA continuum is assumed to carry 
points of recombination, Y, which may result from a fairly compli- 
cated scheme such as (1) independent replication of separate parts of 
the continuum followed by synthesis of these parts into one linear 
whole, or (2) a sequence of “copy-choice” alternatives between two 
parental “continua” or (3) both (1) and (2). The “genetic distance” 
between loci in a virus would then be the average number of points 
of recombination, Y, carried between the two loci in the virus 
progeny. The frequency of recombination between these loci is 
asymptotic to the probability that an odd number of recombination 
points occurs. However, this is of no value in the determination of the 
true genetic distance unless the distribution of the number of recombi- 
nation-points in this region per virus offspring is known. The theo- 
retical determination of this distribution from a suitable model 
concerning the distribution and number of “copy-choice”’ alternatives 
or recombinant syntheses of replicated parts, which underly the 
occurrence of Y, is complicated and has not yet been worked out 
satisfactorily. 

In the case of classical single-strand analysis, the event Y, is the 
cross-over producing two reciprocal recombination points on homo- 
logous chromatid strands at homologous sites. As before, the proba- 
bility that a strand carries a recombination for two genes, is the 
probability that an odd number of recombination points occur on 
the strand connecting the two genes. Generally speaking, this proba- 
bility is 

oo 
A-B) = #1 — Xi (P — 9)” fan), (1.3) 


n=0 


where fay is the probability that exactly 2n crossovers occur on the 
4-strand set between them, and p and g are the probabilities that 
adjacent crossovers form 2- and 4-strand double exchanges respective- 
ly). (f p =g, (A-B) — (1 — fo)/2, PERKINS (1954)). If values of 
b,q, fo, f2, .…. were known, it would be possible to construct a function, 
g, not necessarily single-valued, such that 


1) It is not necessarily true that p and q remain constant. In fact, the data 


of LINDEGREN and LINDEGREN (1942) indicate that they may vary from region 
to region. 
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À = g(f(A-B)) = distance = Y nfn. (1.4) 
n=0 
This function, g, is called a “mapping function”. A number of em- 
pirical functions have been obtained (notably by HALDANE, 1919, 
and KosAMBr, 1944), and functions for special values of the unknowns, 
for example 


In (1 — 2f(A-B)), (1.5) 
where p — g, fn == Ate-À/n! (TRow-HALDANE) 
and the inverse of 

KA-B) = HI — eÀcosh (AVp — 4) (1.6) 


where p and q are arbitrary, and 


Dn AneAfn! 
(SHULT and LINDEGREN, 1956). 


On the basis of the Trow-Haldane formula, it is easy to show that 
if p —gq, and the sequence (f„) follows a Poisson distribution, the 
coincidence index will be unity 4). 

Cases in which p:g:r #1:1:2and f, # AreÀfn! are termed 
chromatid and chromosomal interference, respectively. The presence 
of either of these conditions usually leads to nonunitary coincidence 
values 2). 

When the coincidence index is unity, the Trow-Haldane formula 
yields additive distances, and interference is normally assumed to be 
absent. When this formula yields distinctly nonadditive distances, 
there are two possibilities: (1) the genes are written in the incorrect 
order, or (2) interference of either type is present. 

Further, there are only two topologically distinct modes of non- 
additivity in all three possible gene orders: 


1) Note that the relative frequency of 3-strand double exchanges, „ has no 
effect upon the coincidence index. Further, the assumption of a Poisson distri- 
bution presumes an infinite set of possible cross-over sites on the continuum. 
This distribution would hold approximately if there were a large but finite 
number of potential sites of crossing-over. 

2) This is not always the case, however. In the very interesting model of 
SANDLER and KASTELBAUM (1958), coincidence is unity even though the distri- 
bution of the number of exchanges is not Poisson. This proves that unitary 
coincidence is not evidence for the absence of chromosomal interference. 
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Case 1: Coincidence is less than or equal to unity in one order but 
decidedly greater than unity in the two remaining orders, or Case II: 
Coincidence is greater than unity in all three orders (see figure 1). 


COINCIDENCE INDEX 


GENE ORDER CASE I CASE II 


RELATIONSHIP BETWEEN DISTANCES 


FROM FROM 
NON-ADDITIVE NON-ADDITIVE 
© 
B 
Ep 
TO td 
NEARLY ADpitive || NEARLY SDE 
ee e—*—6 


THE TWO POSSIBLE DISTANCE RELATIONSHIPS IN A 
THREE-POINT TEST 


FIGURE | 
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Both cases have been observed in Drosophila, although Case I 
predominates. 


The equation: 
(A-C) = (A-B) + B-C) — 2(A-B)f(B-C)cane (17) 


(where can — coincidence index obtained for the assigned order 
A-B-C) holds irrespective of the true gene order. (Basically, the 
equation (1.7) depends only on the definitions of the symbols f, 
caBo, and the fundamental property that two recombinations cancel 
each other). Two other equations involving cBca and ccag can be 
obtained from (1.7) by permuting the symbols 4, B, and C, and these 
equations also hold irrespective of the “correct” gene order. It can 
be shown that the coincidence index cannot exceed unity in more than 
one of these three equations. This proves that cases I and II exhaust 
all topological possibilities. 

In Case I, the Sturtevant criterion states (by definition (b)) that 
the genes are in the linear order 4-B-C (see figure 1) with coincidence 


index less than unity. This view is simpler than the supposition that 
the genes are in the order 4-C-—B with coincidence index greatly 
exceeding unity. Thus, there can be little doubt of the correct gene 
order in Drosoplla, where the gene order is not only confirmed by a 
large number of three-point tests, but many of these tests involve a 
coincidence index less than unity (Case I). 

However, Case II where all coincidence indices exceed unity, 
presents two types of difficulties. First, three genes in the order 
A-B-C may involve a region or combination of regions on the 
chromosome in which (negative) interference is so intense that a 
coincidence index much greater than unity is produced. When the 
correct gene order is not known this situation could lead to the 
assignment of an incorrect order to the genes on the basis of the 
Sturtevant Criterion. Second, if three inheritable, Mendelian-segre- 
gating factors are not in a linear order, the criterion would lead to 
the assignment of a linear gene order, when none existed, and the 
consequent “discovery” of negative chromosomal interference results 
from the observation that the coincidence index is now greater than 
unity. Thus, interference actually becomes a “correction” for the 
nonadditivity present in the assigned gene order !). 


1) For problems involving Case II in virus genetics see EDGAR and STEINBERG 
(1958) and DOERMANN and CHase (1958). 
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3. Criteria of Gene Order in Tetrad Analysis 

Three-point tetrad data contains 12 tetrad classes and possesses 11 
degrees of freedom (when both kinds of 3-strand double exchanges are 
distinguished) thus yielding far more information than either two- 
point tetrad analysis, or three-point single-strand analysis. With the 
exception of three special cases, such an analysis represents the only 
observation system which allows a definitive distinction between 
chromosomal and chromatid interference !). 

As a consequence of the larger amount of information available 
from three-point tetrad analysis, there are two independent methods 
for determining the order of loci. 

Whitehouse Correction. Tetrads from three-point analysis fall into 
five classes: (1) nonrecombinant, (2) single recombinant, (3) double 
recombinant, (There are two types: (a) a single recombination in each 
of two adjacent regions, and (b) a four-strand double recombination 
in one region), (4) triple recombinant, consisting of a 4-strand double 
recombination in one region, and a single in the other, and (5) 
quadruple recombinant, consisting of a 4-strand double exchange in 
each of the two adjacent regions. When regions and strands are 
distinguished these five classes break up into the 12 classes previ- 
ously mentioned. It is clear that in a three-point analysis, the 
recombination classes mentioned above do not represent all ex- 
changes because of the bias introduced by the fact that a 4-strand 
double is the only double exchange detectable within a single region. 

The Whitehouse correction takes this bias into account. Let us 
define the chromatid interference vector as the vector (bp, 7, q) where p, 
r and qg represent the probability that adjacent exchanges (that is, 
those exchanges not separated by intervening exchanges) will be in 
a 2-strand, 3-strand or 4-strand relationship respectively. Let the 
chromatid recombination vector be denoted. (P, R‚Q) where P, R and 
Q are the respective probabilities that double recombinations (formed 
from single recombinations occurring in each of two adjacent regions, 


1) Multipoint analyses of single-strand data fail to be definitive here since 
region-to-region variability in the intensity and type of chromatid inter- 
ference producing region-to-region variability in coincidence values, has never 
been excluded (WEeINSTEiN, 1958). Tests for distinguishing chromosomal and 
chromatid interference in single-strand and tetrad analysis see SHurT and 
LINDEGREN and SHULT (in ms.). 
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class 3-a above) exhibit a 2-strand, 3-strand or 4-strand relationship. 
The Whitehouse correction allows computation of (P,7,q) from the 
twelve 3-point tetrad classes. It is based on the following assumptions: 
(1) (®, 7, q) is constant for all exchanges in the three-point interval, 
and (2) triple exchanges do not occur in any sub-region, nor do 
gquadruple exchanges occur in the total three-point interval. 

It can be shown from Whitehouse’s equations (WHITEHOUSE, 1956), 
and, in fact, more generally for an arbitrary distribution of exchanges 
(SrurT, 1956) that |p —g| > |P —Q| and [Ir —05| > IR — 0.5]. 
These inequalities imply 


Xb, 7,4) > XP, R, 0), (2.1) 


where 


X2(a, b, c) — Ala — 0.25)? + 4(b — 0.25)2 + 2e — 0.5)? =— 
— 4(a? + b2) + 22 +3. 


That is, (p : 7 : q) deviates (by the criterion of total squared differ- 
ences) from 1 : 2: 1 more than P : R:Q does. Assumption (2) of 
the Whitehouse Correction may be checked to some extent by 
examining Class 5 tetrads, the apparent quadruple exchanges. The 
Whitehouse Correction is intended to apply only when this class is 
empty. However, when an excess of 2-strand double exchanges is 
involved in the application of the correction it is possible for this 
class to appear empty even when a number of quadruple exchanges 
are involved. Thus, at best, assumption (2) leads to an underestimate 
of the number of exchanges. Further, if quadruple exchanges are 
involved, the Whitehouse estimate (p,7,q) then deviates even less 
from (4, 3, 4) than the actual values. That is to say, the Whitehouse 
correction is an wnderestimate of the intensity of chromatid inter- 
ference if assumption (2) does not hold. Similarly, if assumption (1) is 
denied (STRICKLAND (1958) holds reservations concerning this as- 
sumption), the vector (p, 7, q) would be heterogeneous over the three- 
point interval. (In fact, the data of LINDEGREN (1937, 1942) and data 
presented in the final section of this paper indicate that inasmuch as 
(P, R‚Q) varies with the regions considered, the actual chromatid 
interference vector probably is heterogeneous in many instances). 
In this case, the Whitehouse estimate would represent a kind of 
average chromatid interference vector for the entire three-point interval. 
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Thus, the true values, on the whole, would have to deviate from 
(E, 4, 4) even more than the Whitehouse estimate. Therefore, if as- 
sumption (1) is denied, the Whitehouse Correction again becomes an 
underestimate of the true chromatid interference vector. This is to 
say, whatever misgivings one may have regarding the assumptions 
underlying the Whitehouse Correction, the true interference vector 
still deviates from random at least as much as the Whitehouse estimate 
does. 

Method I. The equations of the Whitehouse Correction, in addition 
to giving an estimate of chromatid interference, also provide a means 
of estimating the distribution of the number of exchanges in a three- 
point interval. Thus, once (p,7,q) is determined, the number of 
nonexchange, single, double, triple, and quadruple exchanges, #0, #1, 
na, nz and #4, can be determined from the apparent exchange classes, 
1, 2, Sa, 3b, 4 and 5 listed above, by the equations: 


2 
vrtidlgean (22) 
q q° 
edible ACE eest 
q q qe 
na — [32] 4 (: | vinnen 
q q q° 
Ì 2 
Rene tee 
q q q 


edi 


2rq l 2r 
n= (1 + \ia Ana 

q q q 
Ee (1 d ie) 


where [k] is the frequency of class £. 


In this respect the Whitehouse correction also allows correction of 
the distribution of the nwmber of apparent exchanges in classes 1, 2, 
3, 4 and 5 as well as the distribution of kinds of doubles. Thus three- 
point tetrad data can be used to order loci on the criterion that the 
preferred gene order is defined as that order which exhibits the least 
number of multiple recombinations. This procedure is analogous to 
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classical method of arranging loci from three-point single-strand data 
(STURTEVANT, 1916) 1). 

If three loci are artificially arranged in the incorrect order, some 
exchanges will lead to the observation of recombination in two 
apparently different regions. Thus all incorrect orders will lead to 
higher frequencies of multiple recombination. There are several 
other methods, less refined than the use of the Whitehouse correction, 
which produce a corrected distribution of the number of exchanges 
in a three-point interval. Although these methods are more crude, 
they suffice perfectly well in choosing between “correct” and 
“Imcorrect'’ gene orders. One method involves the simple assumption 
that chromatid interference is absent. Then each four-strand double 
exchange occurring within a single region simply represents 4 of all 
double exchanges occurring within this region. One-fourth of these 
would have been mistaken for nonexchange regions, the remaining 
half (3-strand doubles) being previously classified as single-exchange 
regions. By making the necessary corrections for both regions of the 
3-point interval a new distribution of the number of exchanges can 
be obtained 2). The total number of exchanges, 7 = ni + 22 + 
+ In + Iz + 44, may be calculated from this distribution. The 
correct gene order is then that which yields the smallest value of 7. 

Method II. A further test of gene order available from three-point 
tetrad data can be obtained from the observation that any exchange 
producing recombination in two regions defined by improperly ordered 
genes should involve the same strands of both regions. That is to say, 
when genes are written in the incorrect order, they tend to produce an 
excess of apparent 2-strand double recombinations (P > Q), in 
addition to the excess of multiple exchanges noted in the previous 
paragraph. Thus one may use the criterion that the “correct” gene 


1) It should be kept in mind that the basic assumption that the locí in 
question are in some linear order drastically limits the possibilities of ar- 
rangement. Therefore, neither the method used by STURTEVANT nor the 
analogous method stated here can properly be said to prove that genes are 
in a linear order. Rather, these methods provide an operational defimition 
of order. It is the fact that when a large number of loci are analysed in this 
way, no contradictions in order appear, which indicates that the operational 
definition “works”. 

2) This can be accomplished, by putting p = 4,7 = &, and g = }in equations 
(2.2). 
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order is that order which exhibits the least number of apparent 2- 
strand double exchanges. This method is independent of method 1, 
and it is conceivable that both methods could give contrary results. 

It is important to note that observed instances of an excess of 2- 
strand double recombination (e.g. LINDEGREN, 1942; WHITEHOUSE, 
1958, etc.) cannot be due to faulty gene arrangement since they show 
a large excess of 2-strand double recombination in al/ arrangements. 
The chromatid interference reported in each case, though significantly 
different from random, represents the least excess of 2-strand double 
recombination for all gene orders possible. 


Tasre II. Tetvad Distributions of Linkage Group I 
Uhh HL STSEEN 


HYBRID 108 EE 
ur CH ht 1s an 
Tetrad 
nne frequency 
Distri- en 
Markers bution *) N12 Pi Xa? Ps X32 P3 f 
recombi- 
1 Lal mr nation 
AN-—ch 40 16 95 | 10.3 0.6% 12.4 O2 0.96 | “3.50% „4205 


ANSTEIN MO 2E ZA NS oel O A F2O88 zal Ons 468 [<< 1OS6 ‚2314 
ANT She al Sen 0 SSN LE 2E lr 152 Olmiesel0e Shall 


ZUNE ADN SIN ZEN ZS 31.0% 29 14,0% 0.82 | 37.0% „4830 
ch_-HI 62| 6| 84 | 46.1 One 70.8 ie KOS 8.9 0.28% ‚3158 
ch-IS 441 12| 97/ 18.3 O4 20.8 |4 Xx 105 0.74 | 42.0% „3954 


ch-UR SOU SZ) N22 LO ALE OSI LOZIE OS „2455 
HI-IS |126|\ 2| 47| **) [6.45 x 10-35 | 388.4 silO ee hi ZA De „1457 
TELESIS SO 015% 12.3 0.21% De en 4432 
IS-UR Sietze eZ Al 4,3% LAS 26 miel 15% „4598 


*) Occasional tetrads involving irregular segregation of one of the markers are excluded 
from this compilation. 

1? = (IL — ID2/(I + II), Pi = P(X) (laf), 

X2? = value of Chi-square obtained from a test of the hypothesis, 1: II: III = 1: 1:4 
Pa= Plan (2rdet) 

X3* — value of Chi-square obtained from a test of the hypothesis (I + II) IIS 2 : 4 
IER TAN (OL le. 

**) In these cases Pj was calculated on the basis of a binomial distribution since the 
number of type II tetrads was less than five. 
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TABLE III. Relative Error in Additivity of Map Distances 


(Error expressed in per cent) 


Definition of Symbols for Different regions 


ee En GE HI IS AN 
Ri Ra Rs 
Ri2 Raza 
R23 
R123 
Raza 
R1234 


Addition of Regions 


Rie — Ri + Pa 
Raa == Rs + Ra 
R23 — Ra + R3 


Ries — Ri + Re + P3 
Ries — Rie + Pz 
R123 — Ri + Ra3 
Raza —= Ra + R3 + Ra 
Raza — Raz + R4 
Raza —= Ra + Raa 


Ri234a = Ri + Ra + R3 + Ra 
Ri234a = Ri + Ra + Lz4 
R1234 = Ri ze R23 td Ra 
Riza = Ria + Pz + Pa 
Riesa — Ria + Raa 

Ri234 — Ri23 + Lá 

Ri234 — Ri + Roza 


Average absolute relative error 


ke) 

% Error of 

distance in 
centimorgans 


— 26.64 
— 10.97 
— 16.71 
—53.76 
— 28.07 
— 39.38 
—36.17 
— 20.45 
—30.13 
—69.37 
—64.12 
— 55.69 
—44.92 
— 39.66 
— 18.19 
— 37.88 


+37.01 


% Error of 
distance in 
Stranes 


23.18 
3.96 
14.18 

19.91 
—0.16 
14.00 
13.27 
111 
11.93 
32.86 
32.13 
26.30 
17.93 
17.20 
18.02 
25.65 


16.99 


% Error of 
distance 
measured 
by (2.1-—2) 


Zeo 
— 11.07 
— 0.48 
— 0.19 
— 5.01 
+ 0.51 
— 7.63 
— 7.26 
— 3.05 
— 4.81 
— 1.55 
— 4.55 
TO 
— 6.62 
— 4.97 
+ 0.58 


4,74 
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4. Linkage Group I of Saccharomyces cerevisiae 

Two-point tetrad data of hybrid 108 are presented below in 
table II. The map distances in figure 2 were obtained on a basis of 
(a) Centimorgans — 50f, and (b) Stranes, distances obtained from 
the Haldane-Trow formula (D — —50 In (ai — 42), see (1.5) and 
LINDEGREN, SHULT and DESBOROUGH, 1959) and (c) the formulae: 


A= 5In eene) (2.1) 


41 — 42 
and. 


Distance — Ae0-(aj — 42), (2.2) 


based on Perkin’s model of crossing-over (BARRATT, et. al., 1954). 
(See Appendix I for derivation). 

Maximum additivity is attained with the gene order UR-CH-HI- 
IS-AN. Although the additivity is not perfect (the average absolute 
relative error !) is 37.01% (see table III) for Centimorgans, 16.99% 
for Stranes, 4.74% by (2.1) and (2.2)), any other gene order leads to 
disastrous conflicts in additivity by any of these mapping functions. 
Further, if we assume that the “correct” mapping function (whatever 
it is) is a monotonic increasing function of f,, this gene order produces 
no instances in which a subregion has a distance longer than that of 
a region enclosing it. 


5. Application of Ordering Criteria to Three-point Tetrad Data of Linkage 
Group 1. 


Clearly, since hybrid 108 was heterozygous for five loci of linkage 
group L, there are ten different combinations of three genes for which 
three-point tetrads may be collected. These are presented in table IV. 
In these data both kinds of three-strand doubles are not distinguished 


1) Note that total distance wnaderestimates the distance expected by adding 
sub-regions when distances are calculated as Centimorgans, whereas the total 
distance tends to overestimate in some regions hwen Stranes are used. This 
suggests that the correct formula lies somewhere between the two estimates 
and that either intra-regional positive chromosomal interference or inter- 
regional positive chromatid interference (excess of 4-strand doubles) is present. 


Three-point tetrad analyses presented in this next section definitely rule out 
the latter hypothesis. 
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UR-CH-HI UR-HI-AN 
HI-AN 
Led II 
Taalikel 14 
UR-HIS I | 5 | O 2 
mer o [ie Jie Ta 
UR-IS-AN 
 # IL iik ï ree man EN 
I [30 ne 2 I lislo Sen 
CH-HIS {3 |O UR-Is< IH |6 |O 3 
ENEN ie ojee 7 
HI-IS-AN UR-CH-AN 


II I 

1 [21 [o 25 

CH-HI& JI | 2 | O | 
I{sojojs ls |4 


UR-HI 
UR-CH- IGE 
CH-IS bern CH-IS-AN ET 
badk Ruf Ee II 
Ij9 |2 I4}25 | 0 9 
Et Ee lkhed sla APD) e) 
Heal 2 mj4aalo)3s [5 [e 
TABLE IV 


so that only 11 (rather than 12) tetrad classes are distinguished. The 
Roman numerals heading the rows and columns of each box denote 
the tetrad type for the two genes listed over the rows and columns 
respectively. Thus, for the triplet UR-CH-HI, there were 11 tetrads 
which were type I (parental ditype) for UR-CH and CH-HI simul- 
taneously. When a 3-point tetrad is type III (tetratype) for both 
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pairs of loci concomitantly, the two reciprocal recombinations 
occurring in each sub-region (e.g. UR-CH and CH-HI in UR-CH-HI) 
may be arranged as either a 2-, 3-, or 4-strand double recombination. 
For this reason, the “III-III”’ box is divided into three smaller boxes, 
the entries denoting the number of tetrads that are respectively 2-, 3- 
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Fig. 2. Three maps of Chromosome I, obtained by the application of three 
mapping functions (a), (b) and (c) (see text) to data of table LI. 
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and 4-strand double recombinations when read from left to right. 
Thus in the triplet UR-CH-HI the double recombinations, formed 
single recombinations occurring in each sub-region, are comprised of 
five tetrads that are 2-strand doubles, six tetrads that are 3-strand 
doubles, and two tetrads that are 4-strand double recombinations. 
In presenting these data, a specific gene order is presumed for the 
purpose of defining the two sub-regions indicated on the row and 
column margins. In each case the two sub-regions are contiguous (if 
not overlapping) for they contain a common marker, and the order 
is that given in the map in figure 2. However, presentation of the 
data as if the specific gene order given above each set of boxes were 
“correct”’ introduces no bias, since a particular permutation of the 
box-entries would allow the data to be presented equally well if the 
genes were in a different order. 

Tetrads involving irregular segregation of one of the three markers 
are not included in the data of table IV, since they do not admit 
unambiguous classification. Further, tetrads in which only three- 
spores survived are not included in this table even though it would 
be possible to infer the genotype of the missing spore on the assumption 
that only regular segregation is involved. If irregular segregation 
occurred for one marker in a tetrad in which only three-spores 
survived a “recombinant” strand not due to crossing-over would 
result but would be inferred as a reciprocal recombinant, thus 
biassing the tetratype frequency. It would not be valid to use such 
a tetrad in a critical analysis of gene order since recombinations of 
this kind (perhaps resulting from gene conversion, post meiotic 
mutation, or possible segregation of a particle-factor) do not neces- 
sarily involve actual recombination points on the ordered continuum 
of the type “Y”’ which are fundamental to the principle of order (see 
discussion in section 2). Thus table IV contains only those tetrads in 
which all four spores survived and revealed regular Mendelian 
segregation for all three phenotypes. For this reason, the marginal 
frequencies of each diagram in table IV will not coincide with the 
two-point data (tetrad distributions) of table Ill, nor coincide 
exactly with marginal frequencies for the same region in other 
diagrams of table V. 
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Fig. 3. Distribution of the number of multiple exchanges for three possible 
gene-orders of three different gene triplets. The correct order in each case 
(top row), shows the least number of multiple exchanges, and fits the Poisson 
expectation more closely. The preferred gene order also shows the least number 
of 2-strand double exchanges. The ratio above each histogran denotes the 
count of 2-strand, 3-strand, and 4-strand double recombinations for that 
gene order. Multiple exchange distributions were corrected in order to account 
for double exchanges occurring within a single region. This correction is com- 
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TABLE V. Gene Orders Uniquely Determined by Method III 


(1) UR-CH-HI 

(2) CH-HI-IS 

(3) HI-IS-AN 

(4) UR-CH-———-IS , confirming (1) and (2) 

(5) CH-HI___—AN confirming (2) and (3) 

(6) CH———-IS-AN confirming (2) and (3) 

(7) UR-CH-———_———— AN confirming (1), (2), and (3) 


6. Application of Method I 

The three-point data of table IV may be used to order these loci by 
method 1. For illustration consider the triplet UR-CH-HI. There are 
ll nonexchange tetrads, 29 + 33 —= 62 single exchange tetrads, 
142564 2== 16 double exchange tetrads, 1 + 2 == 3 triple 
exchange tetrads and no quadruple exchange tetrads; these are the 
“apparent” exchanges within UR-CH-HI. However, three of the 
doubles are 4-strand double recombinations within a single region, 
indicating that three of the nonexchange tetrads and six of the 
single-exchange tetrads are actually 2- and 3-strand doubles, re- 
spectively. Similarly, the three tetrads which are apparent triples 
each involve a four-strand double recombination within a single 
region, indicating that three of the singles and six of the doubles should 
also be classified as triples. These corrections are diagrammed below: 


No. of exchanges e) 2 ge) 4 


Distribution of 
“apparent” exchanges 11 62 16 3 0 


Corrections —3 | —6 | +9 — 
RÓ  Ö 
Corrected distribution 8 53 19 12 ©) 


of exchanges 


puted from the observed frequencies of 4-strand double exchanges within a 
single region and the assumption of no chromatid interference. This correction 
may be refined by supposing the chromatid interference to be the same as the 
double recombination ratio given above each histogram. When corrected in 
this way, some gene arrangements lead to negative multiple exchange frequen- 
cies (usually in the no. of singles). These arrangements are marked by an 


asterisk. 
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The total number of exchanges involved in the corrected distri- 
butionis 53 E21 ZZ 

Figure 3 shows a series of histograms denoting the corrected number 
of non-, single-, double-, triple-, and quadruple-exchange tetrads for 
a triplet of three genes arranged in each of three possible orders. Note 
that in each case, when the triplet is written in the incorrect order, 
the total number of exchanges, 7, increases drastically, due to a rise 
in the frequency of multiple exchange. In particular, note that those 


classes involving an even number of exchanges (zero, two or four) 
are favored beyond the Poisson expectation (LS-CH-HI, [S-UR-HI, 
AN-UR-HI, AN-UR-IS and AN-CH-IS are extreme examples). 
This is precisely what would be expected if one part of the chromosomal 
region were being counted twice because of the incorrect ordering of 
the loci. Note that by the Sturtevant Criterion (minimum number of 
multiple exchanges) the same order obtained by examination of the 
2-point tetrad distributions (in the previous section) is confirmed. 


7. The Application of Method 1I 

The second criterion (Method II), which defines the “correct” or 
desirable order as that arrangement of genes which leads to the 
smallest number of apparent two-strand double recombinations 
confirms the Sturtevant Criterion in each case. Thus, with the three 
orders, UR-CH-HI, CH-HI-UR and HI-UR-CH, the chromatid 
recombination ratios (2- : 3- : 4-strand double recombinations equal 
P:Q:R)are5:6:2,29:6: 2, and 33: 6: 1, respectively. Since 
the first order involves only five 2-strand double recombinations, 
UR-CH-HI is regarded as the “correct” gene order for these three 
genes. It can be seen in Figure 3 that the minimum number of 2- 
strand double recombinations parallels in each of the ten cases the 
minimum value of 7, the corrected total number of exchanges. 


8. Method III: A Refinement of Method I 

It is possible to refine Method I by replacing the assumption 
(bp, 7, 7) equals (4, 5, &), by (p, 7, q) equals (P, R‚ Q) in equations (2.2). 
Thus, in the case of the gene-order HI-IS-AN, a four-strand double 
exchange represents two fifths, (2/(O + 3 + 2)), of all double ex- 
changes in a single region. This will lead to a slightly different cor- 
rected distribution of multiple exchanges. Thus from table V the 
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distribution from tetrads among the apparent exchanges classes is: 
69 nonrecombination tetrads, 46 singles, 5 doubles, 1 triple and no 
quadruples. The single triple-recombination tetrad involves a four- 
strand double recombination in one region (HI-IS). Since for each 
four-strand double there are at least 1.5 3-strand doubles, but no 
2-strand doubles, the corrected distribution becomes 69-46-3.5-4.5-0. 
This method yields a slight refinement of the distribution 69-45-3-5-0 
obtained by method II. However, when this refinement is applied to 
the incorrect order, extreme frequencies of multiple exchanges result — 
sometimes even leading to negative frequencies for some exchange 
classes. The gene order /S-AN-HI (which involves the overlapped 
interval /S-AN counted twice), there are eighteen apparent 2-strand 
doubles, three 3-strand doubles, and only one apparent 4-strand 
double. Thus each four-strand double within a single region represents 
only one out of 22 double exchanges, making the corrected distri- 
bution of multiples exchanges, 69 nonexchange tetrads, -8 singles, 
16 doubles, 44 triples and no quadruples. The gene order AN-HI-IS 
also leads to negative frequencies when corrected in this way: 69 
nonexchange, -43 single, 25 doubles, 68 triples, and no quadruples. 
Some gene orders yield frequencies that can be described as nothing 
short of ridiculous. For example, the corrected distribution for the 
gene order CH-HI-IS is 25.5-38-18.5-9-0, while the distribution for 
HI-IS-CH is (-6)—(-10)-56-51-0, and that for 1S-CH-HI is even 
(-3,645)—(-930)—(-230)-94-4,797. The latter two distributions imply 
that for these gene orders, setting (bp, 7, 7) = (P, R‚ Q) is an incorrect 
assumption; that, in fact, X2(p, 7, q) < X?(P, R‚ Q), contrary to (2.1). 
If these were correct gene orders we would have to suppose that 
adjacent exchanges close together in the same region do not “inter- 
fere”’ (in the sense of chromatid interference) as much as exchanges 
further apart, on the average, in different regions. 

The logic behind this is as follows: Consider the assumption (e) 
that chromatid interference is effected only through adjacent exchanges. 
We may suppose also (f}) that the three genes under consideration 
exist in some linear order. The assumption (@) then takes three forms 
(Bi), (B2) and (B3) for the three possible orders of these genes. The 
assumptions («) and (@;) taken together imply that method II is an 
under-correction, i.e., quantities subtracted from class frequencies in 
the process of correcting are less than the “true” amounts to be 
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subtracted. Thus the assumptions imply that negative frequencies 
cannot occur in the corrected distribution, and we write: 


(a) & (Bs) 9 nonnegative frequencies (4 — 1, 2, 3). 


Denial of the right hand side of this implication forms the contra- 
positive, and implies not-{[(«) & (p‚)], that is, either not-(«) or 
not—(f). If we take assumption («) as true, then the existence of 
negative class frequencies resulting from this method of correcting 
the distribution of apparent exchange classes for a particular gene 
order, is demonstration of the impossibility of that gene order. In 
particular, the gene orders AN-HI-IS, HI-IS-CH and IS-CH-HI 
mentioned above, must be impossible since they were associated with 
negative class frequencies. All gene orders impossible by this criterion 
are marked with an asterisk in Figure 2. When only one order is possible, 
this is regarded as the “correct” order defined by this criterion. 

In three out of ten instances, the test for order by Method [II 
is not clear-cut. In one instance (UR-IS-AN), the method does 
not yield a wmigue “correct” order, and in two cases (UR-HI-IS 
and UR-HI-AN) no arrangement is possible. Here we have 
not-[(a) & (B1)] & not-[(a) & (B2)] & not-{(a) & (B3)] which implies 
not—(a) or not—(f) or both. Thus the above two cases are tantamount 
to the assertion that either assumption («) is false or the genes are 
not in a linear order at all. (Both assumptions will be shown to possess 
exceptions in later sections of this paper). 

The remaining seven cases yield unique correct orders. The first 
three of these uniquely define the total arrangement UR-CH-HI-IS- 
AN((1), (2), and (3) in table V). The last four confirm the first three 
in that these orders are deducible from the first three. 

Even though two of the possible gene orders in each 3-point test 
would be rejected either by the multiple exchange criterion of 
method IL, or by the 2-strand double exchange criterion of Method II, 
some of these gene orders are impossible for reasons other than either 
of these criteria, namely by the appearance of negative exchange 
classes when corrected by a method III. 


9. Zmplications of the Inconsistent Order UR-AN-IS 
Although methods 1 and II agree throughout, the triplet UR, IS 
and AN presents two orders which are almost indecisive by the 
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Sturtevant Criterion (UR-IS-AN, T equals 177, and UR-AN-IS, 
T equals 166). It is, in fact, this triplet only, which yields an order 
inconsistent with the remaining nine tests. If hybrid 108 were heter- 
ozygous for these three loci only, the incorrect gene order UR-AN-IS 
would be assigned on the basis of methods of 1 and II. This shows 
that the application of an ordering principle in any single instance 
does not constitute a proof of gene order in any sense — but rather it 
is a matter of synthesis of a linear gene-sequence from a series of 
mutually consistent gene orders. 

Thus, of the ten tests for order by Methods I and II, three tests 
establish the total order, six tests confirm this order, and one test 
weakly contradicts it. 

Further, there is no way to reconcile the contradiction in terms of 
non-linear arrangements, for UR-CH-HI-IS and UR-CH-HI-AN 
both comprise a linear order without contradictions by these criteria, 
so that the only possible departure from linearity is a branching 
point to the right of H/. But this should appear more obvious with 
the triplet (HZ, IS, AN) than with (UR, IS, AN), contrary to fact. 

Finally, since the contradiction does not result from a non-linear 
arrangement of genes, it must then be due to the manner in which 
crossing-over occurs. This would involve interference (chromosomal 
and/or chromatid interference), which is the real implication of this 
inconsistency. An analysis of interference appears in the following 
two sections. 


10. Test for Intra-regional Interference in Linkage Group I 

If both chromatid and chromosomal interference are absent, 
47 — 42 = (1 — 3a3/2)3/2 so that the expected frequency of non- 
parental ditype tetrads on this hypothesis is given by E(a2) = 
—= HI — ag — (1 — 3a3/2)3/2). (This is substantially the Papazian 
(1951) formula without the approximation). If a2 < £ (aa), then 
either positive chromosomal interference or negative chromatid 
interference may be inferred to be present. (An excess of 3-strand 
double exchanges would also give this result). If a2 > E(a2), then 
either negative chromosomal interference or chromatid interference 
(expressed as an excess of 4-strand doubles or a deficiency of 3-strand 
doubles) may be inferred. (There are further special cases when az 
greatly exceeds E(az), which cannot be attributed to an excess of 4- 
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strand doubles (SruLT and LINDEGREN, 1959). However, these cases 
do not appear in these data and need not concern us here). 
Unfortunately, the calculation of E(az) can be carried out, only if 
a3 < 3. Seven sets of two-point tetrad data, with this property can be 
extracted from linkage group 1; these are presented in table VI. 


Tare VI. Test for Intvaregional Interference by the Papazian Formula 


Pair of Total No. | Observed No. | Expected No. 
Markers of Tetrads of NPD of NPD 4 2 
UR-CH 165 1 7.384 2.4038 1,54% 
CH-HI 152 6 5.323 0.2987 | 77.18% 
CH-IS 153 12 17.750 1.4516 | 14.70% 
CH-AN 151 16 17.259 0.3220 | 74.90% 
HI-IS 175 2 2.051 0.0358 | 97.20% 
HI-AN 175 2 5.726 158S 2 mlniee 

IS-AN 171 0 0.123 0.3507 | 72.64% 


The region UR-CH, which spans the centromere of Chromosome l, 
yields a slight indication of either positive chromosomal interference 
or an excess of either 2- or 3-strand double exchanges in the population 
of cross-overs inside this region (see sec. 3n. l and Appendix II). The 
results are subsignificant, however, and there are no similar indications 
to be found in the other six regions. 

It will be noted that all seven distributions deviate in the same 
direction from expectation — that is, in each case there are fewer non- 
parental ditype tetrads than would be expected on a hypothesis of 
no interference. The total number of non-parental ditype tetrads 
from the seven experiments is 39. The total of the expected numbers 
of NPD tetrads is 55.616, and the probability that a deviation, this 
great or greater would occur by chance in 1142 tetrads is 2.26% per 
cent, which is not quite significant. 

It is possible to devise a test for intraregional interference for the 
three cases in which the tetratype frequency, ag exceeds 66.67%. 
First, supposing interference to be absent, a3 should never exceed 
66.67%. If ag does significantly exceed 3, interference is indicated. 
Table VII shows individual tests, with negative results. However, 
the total of these three experiments yields a significant result: 
381 tetratype tetrads observed, 350.57 expected (two-thirds of the 
total number), z — 2.806 and P — 0.5%. This figure is not reliable, 
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however, since the distributions are not independent. If these five- 
point tetrads were accidentally sampled so that a3 for UR-HI was 
15% instead of 667% (the probability of this, as we have seen, is 
1.92%), then we should expect UR-IS and UR-AN to have tetratype 
frequencies exceeding 3 alsó. That is to say, ag for UR-IS is 72.4%, 
because az for UR-HI is 75%. The results are not independent. 


TABLE VII 
Pair of Observed No. Expected No. 
Markers Tetratype (Elas] = 2/3) i a 
UR-HI 132 LZ SS 2.347 192% 
UR-IS 126 116.00 1.608 10.74% 
UR-AN 123 lS 0.907 36.28 % 


However, this puts us in position to apply a final test for intra- 
regional interference in these three regions. Assuming the more-than- 
two-thirds tetratype frequency to result from sampling deviation, 
the best theoretically permussable estimate of az is exactly 3. Thus we 
should expect aj — az for these three regions. (See table VIII). 


TABLE VIII 
Pair of | No. type I No. type II X12*) P 
Markers tetrads tetrads 
UR-HI 32 12 9.09 02505 
UR-IS 31 17 4.08 4,40% 
UR-AN 32 21 2.28 16.00% 


°) Xi? = (1 — IDT + II). 


The region UR-HI deviated markedly from this expectation, 
demonstrating either positive chromosomal interference or an excess 
of 2- or 3-strand double exchanges (chromatid interference) for this 
region. An analysis of three-point tetrad data for chromatid or 
chromosomal interference, which appears in the following section, 
show conclusively that the significant results obtained for the region 
UR-HI are not due to excesses of either 2- or 3-strand double ex- 
changes but are the direct result of positive chromosomal interference. 


11. Test for Inter-regional Interference in Linkage Group I 
The results of the ten possible three point analyses have already 
been presented in table IV. Chromatid interference is indicated by a 
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significant departure from a 1 : 2: 1 ratio in the observed number 
of 2-, 3- and 4-strand double recombinations. Chromosomal inter- 
ference is indicated, if the various tetrad types are not formed 
independently in each of the two sub-regions of the three-point test. 
A theoretical distribution of the nine combinations of tetrad types 
is computed from the product of the corresponding marginal frequen- 
cies. However, such a distribution does not lend itself to testing by 
Chi-square since nearly every entry involving non-parental ditype 
tetrads has an expected number less than five. For this reason, some 
categories must be lumped. Thus one could reduce each square of 
nine entries into a 2 x 2 Latin Square by subhuming types I and II 
in a common category. The resultant Latin Square possesses but one 
degree of freedom and therefore the test is identical to the usual 
normal-curve z-test of the expected and observed numbers of 
“I-III 3-point tetrads — i.e, all tetrads involving a single reciprocal 
recombination in each of the two sub-regions. A significant deviation 
of the observed number of double recombinations from the expected 
number indicates chromosomal interference and is sharply distinguish- 
ed from chromatid interference in this case. 

The analysis of both chromatid and chromosomal interference are 
presented in fig. 4. Note that there is a general excess of two-strand 
double recombinations throughout the data. The total distribution 
of strand-relationships, 69 : 77 : 31, differs enormously from that 
expected on the hypothesis that recombinations occur independently 
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with respect to the choice of strands involved. However, the signifi- 
cance of this finding is subject to the criticism that the data com- 
prising the total are not independent since common regions are 
involved in these tests. In fact, when only adjacent pairs of simple 
regions are considered no cltromatid interference whatever is detect- 
able. Thus, the total distribution for UR-CH-HI, CH-HI-IS 
and HI-1S-AN, is only 8-13-6, in accord with the expected 1 : 2: 1 
ratio. Although most of the individual three-point tests show sub- 
significant chromatid recombination ratios, the two long regions, 
UR-HI and HI-AN exhibit sixteen 2-strand and only four 4-strand 
double recombinations. This result is significant, with P — 6.196 X 
Xx 2-20 < 0.605%. (Computed from a Binomial Distribution of 2- 
and 4-strand doubles). 

This marks an unusual instance in which exchanges occurring in the 
long region UR-HI tend to be arranged as 2-strand doubles with ex- 
changes in the long region H1-AN, while double exchanges within these 
regions show no chromatid interference. These data contradsct the hypo- 
thesis (a) that chromatid interference is effected through adjacent exchanges. 

It will be recalled that the data of table VIII indicated either 
positive chromosomal interference, or chromatid interference ex- 
pressed as an excess of 3-strand double exchanges and possibly an 
excess of 2-strand double exchanges. Figure 4 shows no evidence of 
an excess of 3-strand double exchanges for the regions UR-CH-HI, 
indicating that positive chromosomal interference must have been 
the major agent contributing to the excess of type III tetrads in the 
region UR-HI. This is born out, somewhat, by the near significant 
observed number of double recombinations for UR-CH-HI: observed 
no., 13; expected number 22.48; z — 2.30; P — 2.14%. Much of this 
effect may be present in UR-CH alone (P = 1.54%, see Table VII). 
The centromere of Chromosome l is located in the region UR-CH. 
Thus positive chromosomal interference is operative in regions 
placed nearly symmetrically about the centromere !). Other regions 


1) It is of some interest to know whether chromosomal interference is actu- 
ally effective across the centromere — i.e. do cross-overs, occurring on one side 
of the centromere inhibit exchanges occurring on the other side? An analysis 
presented in Appendix II indicates that complete interference occurring between 
CH and the centromere is sufficient to account for the observed distributions, 
while complete interference between UR and the centromere is not sufficient. 
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in the right arm of chromosome I show no chromosomal interference 
whatever. 

Taken together, all of these data strongly support the contention 
that interference effects are specific for the combination of regions 
involved, and that observation of chromatid and chromosomal inter- 
ference may be easily obscured in lumped data taken from several 
regions as in Perkins’ (1954) summary of interference data. 


12. A Non-Linear Gene Arrangement 

In family 108, the gene #4 (thiamineless) exhibits less than 50% 
recombination with the genes of linkage group 1 considered in the 
previous section. The tetrad distributions are listed in table LX. 


TaBLE IX. Tetrad distributions of TH and the genes of linkage group 1 


fe Tetrad 
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TH-HI |49| 6| 73 |33.6| < 105 48.7 |< 10-6 5.4 21% 
TH-IS |39| 9| 80 | 18.8 10-5 22.1|2 x 10-5| 1.0| 300% 
TH-AN | 28|13[ 88 | 5.5 2.0% 54A| 70% | O.1| 800% 


*) Pis determined from a cumulative binomial distribution. 


The numbers, X12, X2? and X3? denote the value of Chi-square for 
the hypotheses) L II Solesolpl adlealesdignden 4oand (Bedenk 
: HI == 1 : 2, respectively. An asterix indicates use of the binomial 
distribution when II < 5..Since TH shows its strongest linkage with 
CH, which is internal to linkage group 1, these data strongly suggest 
a non-linear relationship between TH and the genes of linkage 
group [. 

The sequence TH-CH-HI-IS-AN forms a linear map with distance 
measured in Stranes (Dj — —50 In (ai — a?), see fig. 5). This map 
exhibits no contradictions in that no sub-region has a distance longer 
than a region containing it. Additivity for this gene sequence is fair, 
the relative error being 226% (see table X). Thus the sequence 
TH-CH-HI-IS-AN, does not appear to depart from linearity. The 
linear sequence UR-CH-HI-IS-AN was established by the ordering 

“principles of sections 4 through 8. These two sequences can be 
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synthesized into at least one of two possible sequences: UR-TH-CH- 
HI-IS-AN and TH-UR-CH-HI-IS-AN. That is, if linearity holds, 
one of the orders UR-TH-CH or TH-UR-CH must hold. But either 
of these two orders implies a contradiction in distance, for UR-TH 
(which is a simple rather than a compound region in either of these 
arrangements) has a distance of 40.04 Stranes (calculated as 


TH CH HI 


28.19 4 49.90 5 17.24 + er 


| | | 
| 78.23 I | 

| | 
LL | 54.43 stoan) 

| | 
| 91.94 
l | Ì 
| | 72.62 | 
iS EN Tr A | I 
a 107.71 | 


FIGURE 5 LINEAR ARRANGEMENT OF TH WITH GENES 
OF CHROMOSOME I 


TABLE X 
é 9% error of distance 
Addition of Regions red 

Ria —= Ri + Pa ERS 6 
R3a — R3 + Ra — 9.89 
Raz == Rat R3 EED 
Riaz —= Ri + Ra + Ps =D 
Ri23 — Riz2 + Pz + 5.09 
Ri23 — R1 + Ra3 si dE 
Raza — Ra + R3 + La — 26.56 
Raza — Ro3 + Ra — 16.97 
Raza — Ra + R34 220.15 
R1i234 = Ri + Ra + R3 + Pa + 26.05 
Riesa = Ri + Ra + Raa + 29.91 
Ri234 = Ri + Raz + Ra 36.61 
Ri23a — Ri2 + R3 + Ra 34.28 
Ri234 = Ria + Raa + 28.63 
Ri2sa = R123 + Ra +-38.60 
R1234 — R1 + Roza RS 0 


Average absolute relative error 22.6% 
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Di — —50 In (ai — a2)), while the distances for CH-TH and UR-CH 
are only 28.19 and 24.55 Stranes respectively. A similar argument 
holds when distance is calculated from the frequency of tetratype 
tetrads as Dz = —33.33 In (1 — 34a3/2): thus UR-TH (which is a 
simple region in orders UR-TH-CH and TH-UR-CH) is longest, 
measuring 47.57 Stranes while TH-CH and UR-CH are 34.62 and 
42.20 Stranes respectively. The non-additivity in these distances is 
shown in fig. 6. 


40.04 à 47.57 


(a) (b) 
Dr —50 In (a a,) D; — 33.33 In (1-3a,/2) 


Fig. 6. Non-additive distances in Stranes for three genes UR, CH and TH. 

In (a) the distances are calculated from aj — as, the difference between the 

frequencies of parental and non-parental ditype tetrads; while in (b) the 
distances are calculated from az, the frequency of tetratype tetrads. 


It is necessary to test this departure from linearity for statistical 
significance. Since most statistical tests assume normal or binomial 
distribution of the random variables involved, these would not be 
appropriate if applied to the distances involved. For example, 
although the relative frequencies of type IT and type II tetrads (a1 
and a2) may be assumed to follow a normal distribution, the distance, 
measured a —S50 In (aj — az) will not be normal. This difficulty may 
be eliminated by designing the test so that the object tested is a 
tetrad distribution departing at random from a theoretical tetrad 
distribution. 

Consider the three genes UR, TH and CH. If they are in a linear 
order, that order must be one of the following: (a) UR-TH-CH, 
(b) CH-UR-TH or (c) TH-CH-UR. We have already seen, from a 
consideration of the linkage group as a whole, that only (a) and (b) 
are permissable. Let us suppose that (a) is the correct order — i.e. 
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UR-TH-CH. In the absence of inter-regional interference, the 
distances should be additive, so that the expected value of D for 
UR-CH is 40.04 + 28.19 — 68.23 Stranes. Hence the expected value 
of aj — az for the tetrad distribution for UR-CH is E(a1 — 42) — 
— exp (—68.23/50) — 25.9%. Similarly, the expected value of De 
between UR and CH is 47.57 + 34.62 — 82.19. On this basis, the 
expected frequency of tetratype tetrads for UR and CH is E(a3) — 
= 2[Ì — exp (—3D2/2)]/3 = 61.0%. Thus, assuming no interference 
and the linear order UR-TH-CH, the expected tetrad distribution 
for the 165 tetrads of Family 108, segregating regularly for UR and 
CH, is 53.68-10.91-100.41, while the observed distribution was 
85-179. The observed distribution shows stronger linkage than 
could reasonably be expected (X2 —= 32.01, P = ca. 10-7). 

The order CH-UR-TH may be subjected to the same test. If 
distances are additive, the expected tetrad distribution for CH-TH 
is 47.07-7.44-89.48 while the observed is 82-0-62. Again, the observed 
linkage is too strong (X2 — 41.80, P < 10-7). 

In addition to the fact that order (c), UR-CH-TH is inconsistent 
with any synthesis of the linear order UR-CH-HI-IS-AN with the 
additive order TH-CH-HI-IS-AN (the latter is supported by 7 
independent checks of additive distance), the arrangement is also 
non-additive (X? — 15.452, P = 4.5 x 10%). 

Thus all three orders are significantly non-additive. This implies one 
of two situations: (1) The genes are not physically situated in a linear 
order, or (2) negative chromatid and perhaps negative chromosomal 
interference (to depress the tetratype frequency for the distal loci) is 
operative throughout some linear 3-point interval. However, on 
hypothesis (2), the interference must be assumed to be so intense 
that it actually causes the distance between distal loci to appear 
shorter than that of at least one of its subregions !). Negative inter- 
ference of such high intensity has never been reported, even in the 


1) This is required for all possible orders of UR, CH and TH. It is obvious 
in orders (a) and (b). With order (c), UR-CH-TH, and the linear order 
UR-CH-HI-IS-AN established in sections 4-8, the only linear orders which 
can be synthesized are UR-CH-TH-HI-IS-AN, UR-CH-HI-TH-IS-AN, 
UR-CH-HI-IS-TH-AN, and UR-CH-HI-IS-AN-TH, all of which involve 
the requirement that a region appear shorter than at least one of the subregions 


contained in it. 
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small regions encountered in phage (CHASE and DOERMAN, 1958, 
EDGAR and STEINBERG, 1958) or Aspergillus (PRITCHARD, 1955) where 
negative interference appears to be the rule. Thus we are forced to 
the conclusion (1), that the linkage relations of th with other members 


of linkage group l are nonlinear. 


13. Affinity 

Non-linear linkage relations not due to interference, can be recon- 
ciled with what is known of chromosome structure and the process 
of meiosis, by the hypothesis of affinity (WALLACE, 1954, MICHIE, 
1954). It is sufficient to suppose that the centromere of the TH- 
chromosome segregates preferentially with the centromere of 
Chromosome 1, in parental combination at the first division in 
nearly 100% of the meioses that make up Family 108. This prefer- 
ential attraction between centromeres of non-homologous chromo- 
somes at MI has been called “affinity”’. 

If the sites of affinity did not preferentially segregate in 100% of 
the meioses, the NPD frequency for TH-CH would be greater than 
that predicted on the basis of the Papazian formula. But it is less 
than this value (see table 9), so that by Appendix III the frequency 
of preferential segregation is at least 100 per cent. 

The site of affinity, which produces the disparity between PD and 
NPD frequencies of TH and the genes of chromosome I (table 9) can 
be located on a map on the basis of these frequencies. This is ac- 
complished by forming a branched map with distances obtained by 
Di = —50 In (a1 — 42) (compare with 1.5) between TH-CH, TH-UR 
and CH-UR. The branch point of branching marks the site of 
affinity. On this basis it can be seen that there is 2.87 times as much 
crossing-over between UR- and the site of affinity as there is between 
CH and the site of affinity. 

The position of the centromere can be computed by an independent 
method involving the frequency of second division segregation be- 
tween three centromere-markers (WHITEHOUSE, 1957). The frequen- 
cies of second division segregation can be converted to map distance 
by Dz = —33.3 In (l — 3a3/2). The resulting map is presented in 
figure 7 with the centromere of chromosome 1 marked by a black 
circle. The expected position of the site of affinity on this map can 
be determined by remembering that the distance between UR and 
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the site of affinity should be 2.87 times the distance between CH 
and the site of affinity. The resultant position of the site of affinity 
is marked by a black diamond in figure 7 and compares favorably 
with the centromere. 


CENTROMERE 
UR UE 
27.6 3.7 Ios GH 


SITE OF AFFINITY 


TH 


20.0 


DISTANCE CALCULATED BY D, 
FIGURE 7 


In section 12, it was found that distances involving TH were sig- 
nificantly non-additive. As was remarked near the end of section 2 
(appropos Case II), there were two hypotheses to explain this: (1) 
non-linear order and (2) negative interference. Hypothesis (2) was 
rejected on the grounds that the required negative interference must 
make distances of compound regions appear shorter than distances 
of sub-regions contained in them. Hypothesis, (1) on the other hand, 
fits well with what is known concerning “affinity”. This hypothesis 
now gains further credence from an examination of the centromere 
markers « (mating-type locus) and ga-l (galactose fermentation). 
These loci are well established markers linked to the centromeres of 
different chromosomes (LINDEGREN and LINDEGREN, 1949; LEUPOLD, 
1956; SuuLT and LINDEGREN, 1957; HAWTHORNE, 1956; ROMAN, 
1956; DESBOROUGH and LINDEGREN, 1959). Table XI shows F- 
distributions typical of these loci. However in family 108, the tetrad 
distribution for « and ga-l was 52: 27:67. These data show a 
significant excess of parental ditype tetrads over non-parental ditype 
tetrads (X2=7.91, P —0.5%), that is, a preferential segregation 
of centromeres of the « and ga-l chromosomes in parental combi- 


Genetica 
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nation at Meiosis I. This shows that affinity occurs in family 108 
irrespective of hypothesis (1) concerning 7H. The frequency of 


TABLE XI 
Tetrad Distribution for « and gal in eight families of the Carbondale breeding stock 


Family Tetrad Distribution 
93 24: 30: 44 
88 ZORM Nn 
87 21 BLG sn 
80 bles 1 6 

111 Dini 28 SS 
116 39: 44: 57 
90 SR REL ZE 
118 Che seh ed) 


total 153 : 172 : 244 


parental segregation, p, can be obtained from the tetrad distribution 
for « and gal by the formula 


p= HL + ear — 42)] (3.1) 


where Â = —& In (l — 3a3/2) (see Appendix III). The centromeres of 
the « and ga-l chromosomes segregate in parental combination in 
68.5% of all meioses of hybrid 108. 


14. Summary and Discussion 

|. In three-point single-strand analysis, it has been customary to 
assign gene-order on the criterion that for three genes, the “correct” 
gene order is that arrangement which minimizes the double exchange 
class. Thus the frequency of double exchanges being lower than the 
frequency of singles is not evidence of linearity, but rather, the direct 
consequence of the assignment of order. Evidence for linear order 
finds its source in the consistency obtained in a series of assigned 
linear orders involving-common genes. Thus with 4 genes, the two 
orders A-B-C and B-C-D imply the orders A-B-D and A-C-D, 
so that there are two checks of linearity. With » genes, there are 

nl 

(n — 3)!3! 
linear structure increases rapidly as the number of genes are increased. 


— (n — 2) checks of linearity, so that the evidence for 
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2. The greater amount of information available from tetrad data 
allows the construction of three independent criteria of gene order 
(see section 3). 

3. Hybrid 108, which is heterozygous for five loci of linkage 
group Ì was constructed and 186 tetrads were dissected by Mrs. 
LINDEGREN. The genes involved are: UR, uracil-less, CH, cholineless, 
HI, histidineless, ZS, isoleucine-less, AN, requirement for anthranilic 
acid. These five loci provide 10 possible 3-point tests in which criteria 
of gene order may be applied. Thus, with the three criteria mentioned 
above, 30 independent tests of gene order are available. Of these, 
three tests define the order of the five loci as UR-CH-HI-IS-AN; 
22 tests confirm this order; two tests contradict it, and three tests 
are inconclusive. Of particular interest is Method III (section 8), a 
criterion of gene order which presumes that chromatid interference 
(if it is present at all) is expressed entirely through adjacent exchanges. 
For two S-point tests utilizing this criterion, no order was possible. 
This implies either that the genes were not physically situated in a 
linear order (which is unlikely in view of the 22 tests confirming order), 
or that the fundamental assumption underlying this criterion is false: 
that is to say, chromatid interference may operate more intensely 
through non-adjacent exchanges. 

4. This prediction was confirmed in an analysis of multiple re- 
combinations of family 108: For adjacent pairs of simple regions 
ranging from 12 to 50 Stranes in length) there were 8 two-strand 
doubles, 13 three-strand doubles and 6 four-strand doubles, a result 
which does not deviate significantly from a random 1 : 2: 1 ratio. 
However, in three-point tests involving one simple region and one 
compound region composed of two simple regions, the ratio is 
28 : 25 : 13. This result shows a significant excess of two-strand 
double recombinations. This trend is preserved in cases where the 
compound region includes even three simple regions (17 : 23 : 8) and 
in the triplet UR-HI-AN which involves two adjacent compound 
regions (16 : 16 : 4). Thus in these experiments exchanges which are 
further apart or are occasionally separated by intervening exchanges 
actually exhibit a stronger chromatid interference than adjacent 
exchanges. 

5. The total distribution of double recombinations is 69 : 77 : 31 
for 2- : 3- : 4-strand double recombinations, respectively. This fits 
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very closely the ratio 2 : 2: l expected if exchanges are assumed to 
possess parity (i.e. no inate left-right distinction) (Surt, 1957). 

6. Positive chromosomal interference is operative in the interval 
UR-CH. (tables VI, VII, and IX). 

7. The gene f4 (thiamineless) exhibits linkage with the other five 
members of the linear linkage group I, but the result departs sig- 
nificantly from linearity: 


UR 
D, CH-HISTSAN 
TH 


If we define CH-HI-IS-AN as the right arm of Chromosome I, 
then either UR or 7H may constitute the left arm of this chromosome. 
It seems likely that UR is actually the left arm of Chromosome I in 
view of the positive chromosomal interference present between UR 
and CH (otherwise it would be necessary to hypothesize interchromo- 
somal chromosomal interference between the UR and CH arms). 

8. The non-linear linkage relationship may be explained on the 
hypothesis of “affinity”’ proposed by WALLACE (1953) and MrcHre 
(1953). Thus it may be supposed that the centromere of the TH 
chromosome segregates preferentially in parental combination with 
the centromere of Chromosome l at the first division in nearly every 
meiosis of hybrid 108. 

9, The “branch point” of this linkage relationship is shown to 
eoincide fairly well with the centromere of Chromosome I. 

10. As further support for this hypothesis, it is shown that the 
respective centromeres of the « (mating-type) and ga-l (galactose) 
arms segregate preferentially in 68.5% of the meioses of hybrid 108. 

These Criteria of order not only demonstrate the linear order of loci 
of Linkage group 1 for Saccharomyces, but also serve to demonstrate non 
Unearity when this occurs via preferential segregation of centromeres. 


APPENDIX 1 


Derivation of the Mapping Formulae (2.1) and (2.2). 
The distribution of the number of cross-overs occurring per meiosis 


between two genes A and B is given by Perkins’ model as 
eÀ pred 1 (rijn his : 
TE NKT Ee OZR 
St nl St _ nl(ekd— (1 —k) ad 
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where P is the probability that » cross-overs occur. The mean of this 
distribution is 
zn RÀ 


bm iN Pn= É 
et rd EN 


In the absence of chromatid interference, 


k 
Ui me P dn 
1 2 0 TE (3.1) 
and. 
ui — AekMaj — a2), (52) 


which is the formula for distance (average number of cross-overs). 
Putting k — 0.2, (3.1) can be solved for Â as (2.1). (3.2) then yields 
the map distance. (Note that since À is not the average number of cross- 
overs per meiosis, it is not a ‘distance’ parameter, as such. For 
further information on this model, in particular, a mapping function 
for ag based on this model see BARRET ef. al. (1954). 


APPENDIX II 


An Examination of Interference Across the Centromere of Chromosome I 
in Family ro8 

It is known from section 10 that either positive chromosomal inter- 
ference, or an excess of 2- and 3-strand double exchanges occurs within 
the region UR-CH, assumed to lie entirely on chromosome 1. This 
region carries a centromere whose position can be calculated from the 
three tetrad distributions: 


GA-UR 54 53 40 F 
GA-CH 39 31 50 F 
UR-CH 85 1 dj 75 


It would be of considerable interest to determine whether the 
interference operative in lowering the NPD frequency of UR-CH 
actually extends across the centromere. Thus, if positive chromosomal 
interference is operative, do exchanges on one side of the céntromere 
inhibit exchanges occurring on the other side? Or, if chromatid 
interference produces the low NPD frequency, do exchanges on one 
side of the centromere become alligned as 2- and 3-strand doubles 
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with exchanges on the other side of the centromere? This hypothesis 
can be tested equally well by testing its denial. If interference does 
not extend across the centromere, then it must be locallized in one 
or both of the regions on one side of the centromere. If the assumption 
of complete interference in any of these small sub-regions is not 
enough to account for the degree of intra-regional interference in 
UR-CH, evidence that interference extends across the centromere 
results. 

Denial of the hypothesis that interference extends across the 
centromere entails the assertion that exchanges in any one of the 
three sub-regions, UR-centromere, CH-centromere, and GA-centro- 
mere, are independent (with regard both to absolute occurrence and 
the choice of strands involved) of exchanges in any other sub-region. 
On this basis, the frequency of second-division segregation can be 


calculated: 
Marker ed of Second 
Division Segregation 
Ge 7.83% 
On 21.97% 
en 38.31% 


Hypothesis 1 : Complete interference within UR-centromere. No inter- 
ference in CH-centromere and GA-centromere. 

Treating the centromere as a marker, the “tetrad distribution” 
between CH and centromere (where no interference is assumed) is 
0.5912 1 + 0.0257 II + 0.3831 III. With complete intra-regional inter- 
ference in UR-centromere (whether resulting from chromatid or 
chromosomal interference) the tetrad distribution allows no NPD 
tetrads so that the distribution is 0.7803 I + 0.0000 II + 0.2197 III. 

Since our hypothesis entails the assertion of independence between 
sub-regions, we may derive expected values of type 1 and type II 
frequencies for the total region UR-CH and test these with the 
observed. Thus, the expected frequency of type I tetrads for UR-CH 
is 0.7803 x 0.5912 + 0.3831 X 0.2197/4 — 0.4824 and for type II 
tetrads is 0.0257 x 0.7803 + 0.3831 x 0.2197/4 —= 0.0411. Thus we 
must compare the observed distribution of 1 : II = 85 : 1 with the 
expected distribution 79.5 to 6.5. The difference is sub-significant at 
2.56% with X2 — 4,847 (1 df). 
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Further hypotheses, e.g. complete interference in CH-centromere, 
or complete interference in both CH-centromere and UR-centromere 
(though each acts independently) would predict stronger intraregional 
interference in UR-CH and hence be more in accord with the observed 
distributions (P would be well above 2.56%). Even if one were to 
insist that P — 2,56% is significant, the evidence for interference 
across the centromere is not strengthened since interference in CH- 
centromere can easily explain the results. In this case, at most one 
could conclude that interference either extends across the centromere 
or is localized in the region between CH and the centromere. 


APPENDIX III 


Estimation of degree of preferential segregation of centromeres 

Consider two genes situated on different chromosomes. Let us 
assume that the centromeres of these chromosomes segregate in 
parental combination at Meiosis 1 with frequency p and in non- 
parental combination with frequency 1 — p. Let the average number 
of cross-overs occurring between each gene and its centromere be 
Ja and Às, respectively. If interference is absent, and there is no 
preferential segregation at the second division, the tetrad distribution 
for those meioses in which the centromeres segregate in parental 
combination is given by 

Ai — 42 == e-Art)o) 


ag — (1 — e”3rtÂn/2) 


The tetrad distribution for those meioses in which centromeres 
segregate in non-parental combination is the same except that 
a1 — a2 changes sign. By assumption the first case occurs with 
frequency p and the second with frequency 1 — p. 

Thus, the resultant tetrad distribution for all meioses is 


41 — 42 = (2p Ee De Artha) 


while a3 is the same as before. 


Thus,  = (. + Ze Et deel ek 
= 
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INTRODUCTION 


MoRGAN's (1915) famous sentence “The cytoplasm may be ignored 
genetically” probably still represents the attitude of most geneticists 
working with multicellular individuals; nevertheless there are a 
number of papers (there is a good review in MricHarris 1954) which 
claim the existence of “hereditary constituents’’ in the cytoplasm. 
Hence it seems desirable to make a fresh and more systematic ap- 
proach to the problem. Although some authors made a case for an 
entirely nuclear control of quantitative characters in hybrids, it 
appeared to the author that the methods of analysis may have 
obscured possible cytoplasmic effects and that therefore the study 
of hybrids remains probably a sufficiently sensitive measure of these 
possible effects. 

On another side when LERNER (1954) introduced his notion of 
homeostasis, he notes very pertinently; “There are several problems 
connected with an attempt to verify the general idea that hetero- 
zygotes are better canalized than homozygotes. The ideal test would 
consist of comparisons of environmental variability in fitness of 
these types. Unfortunately, practically no data are available on this 
point”. Indeed no systematic study of the influence of environment 
and of its variation on the compared behaviours of homozygotes 
and heterozygotes has ever been attempted. 

Therefore it would seem justifiable: 1°. to devote more attention 
to the general properties of a given nucleo-cytoplasmic combination 
and to become more detached from the rather too scanty phenome- 
nology of hybrids; this is the reason why in our experiments as much 
importance was devoted to the pure parental strains as to the recipro- 
cal hybrids; 2°. and furthermore in order to avoid the danger of 
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considering the nucleo-cytoplasmic relations from a basically wrong 
point of view to make all the crosses in a wide range of environmental 
conditions, i.e. temperatures. 
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MATERIALS AND METHODS 
1. Straus 

A. Diverse laboratory strains have been used: Canton-S, Swedish- 
B, and mutations such as brown (lT, 104.5), radius incompletus 
(III, 47.1), and ebony (III, 70.7) (Bridges 1944). A wild type captured 
in 1957 in Watou (Flandre Occidentale, Belgium) was also used. 

B. +(+), ee(+), ee(e), ++(e) t) (individuals with ++ geno- 
typical formula and (+) cytoplasm, individuals with ee genotypical 
formula and (+) cytoplasm, etc. …). 

The back-crossing scheme shown in table 1 for a Canton-S, + +(+) 
and ebony, ee(e) strain was carried out. After n crosses and two 
generations of selection 4 series are available for comparison: + +(+), 


++H(e), ee(e) and eel +). 


TABLE I 
| (+) cytoplasm strain (e) cytoplasm strain 
P ++) X eed eele)Q Xx ++4d 
B Jel)? and J erossed together | et (e)Q and d erossed together 
F. d(H)Ror Held)L X eed?) | eele)QP) Xx +Hdoretd 
Fz +H(H)P or Held)P X eed?) | eele)Q2) Xx +Hdoretd 


F3, Pa, Fn: the same crosses are mass-repeated for » generations 


1) The two first signs indicate the genotypical formula, the third one between 
brackets, gives the cytoplasm origin. The opportunity and the real significance 
of that special notation will appear in this paper. 

2) Descendant from the Fn-1) generation. 
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C. SS(C)12, CC(S)12 (strain with Canton cytoplasm and Swedish 
nucleus, strain with Swedish cytoplasm and Canton nucleus). 

These back-cross strains were obtained from the Canton-S and 
Swedish-B strains following the crossing scheme shown in table 2. 


TABLE 2 
De 
| (C) cytoplasm strain (S) cytoplasm strain 
P | ce x Ss SQ X Cd 
F, | (C2x SI) x Sq I) (52 x CBP x C4 1) 
Fz | (CQ x SJ)L X SHL Xx Sd) (SR Xx CH) X CH)L Xx CI) 


F3, Fa, Fn: mass-crosses are pursued during n generations 


D. A brother-sister inbred wild type strain was derived from 
Canton-S. 


2. Culture Methods 

A. The Culture Medium. The culture medium is derived from the 
formula: agar-agar 9 grms; cornmeal 60 grms; yeast 6 grms; woodcoal 
3 grms; sugarsyrup 105 cc; water 450 cc; Nipagin N traces. Woodcoal 
is used to give a black color to the culture medium in order to 
facilitate location of the eggs laid. 

The medium was poured into small glass cups 25 mm deep X 26 mm 
in diameter. A glass tube of the same diameter and 85 mm high was 
fixed above the cup and stoppered with cotton-wool. When experi- 
ments involved competition between larvae, the glasscups used were 
of identical diameter and provided a medium surface of 5.47 cm? 
The amount of food per cup was recorded by weighing cups before 
and after filling; only those cups in which the weight of medium fell 
within the range of 11 + 0.15 grms were retained. 

B. Crosses. With a few exceptions, noted below, all crosses were 
made between virgin QQ and 44 aged 4-5 days. Furthermore, these 
flies were derived from the modal range of development time in the 
cultures in which they were grown. In this way a fairly high degree 
of homogeneity among parents would be expected. For each cross 
20-30 pairs were introduced into a single tube for egg deposition. 
Generally enough eggs were laid within a 12 hour period but if this 
was not the case, the flies were allowed to continue laying for a 


1) Derived from the pure laboratory strain. 
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period not exceeding 20 hours. At the end of the 12 or 20 hour period 
the tubes were removed from the cups and the eggs per cup counted 
with a binocular and the number standardized at 300 by removal 
of the excess. 

For each cross three replicates were prepared and transferred at 
17° or 19°, 25° and 30°C for further development. At 17° and 30°C 
the growth was completed in the dark; at 25°C the constant temper- 
ature room was illuminated for 1-10 hours daily. Variations in 
humidity were not taken into account. 

In the denomination of a given cross the maternal parent will 
always preceed the paternal one. 


3. Data recorded 

A. Viability. This was recorded as the number of males, of females 
and the total number hatching from each initial sample of 300 eggs. 
The relatively small surface of medium in relation to the number of 
larvae leads to severe competition. The record viability is thus 
largely a measure of competitive ability. 

B. Development time. The emergence of flies was recorded at 24- 
hourly intervals at 17° or 19°C and at 12-hourly intervals at 25° 
and 30°C. 

C. Development speed. For the given purpose it is desirable to 
calculate the inverse of development time and this is referred to as 
developmental speed or developmental rate (see PRECHT, CHRIS- 
TOPHERSEN and HENSEL, 1955). 

D. Body size. This is based on the thorax length following the 
procedure of ROBERTSON and REEVE (1952), who have noted that this 
dimension is less affected by changes in temperature than other 
possible measures of body size, such as wing length. 

E. Egg production or fecundity (Reirr 1945). Two methods have 
been used: one of them records the life time yield of individual QQ, 
expressed as eggs/day, the other follows the procedure of GOWEN, 
STADLER and JOHNSON (1946) and records the average daily egg 
production over the 6th to the 10th day. When rnales died, they 
were replaced by another one of similar origin and age whenever 
possible. 


1 
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RESULTS 


TI. Relativity of the measured traits as a function 
of the experimental conditions 


In order to avoid as many sources of variation as possible a number 
of preliminary experiments were undertaken to determine an eventual 
influence of experimental conditions such as parental age and parental 
development time, environmental variation, and also to determine 
whether reciprocal crosses differ in relation to such experimental 
variation. 


Ll. Parental age and development time 


Drosopmmla melanogaster from a Canton-S strain have been crossed in 
such a way as to obtain two different crossing groups: in the first the 
development time of the parents is equal in each cross and the parental 
age is different; in the second, parental age is equal but development 
time varies. The crossing scheme is reproduced in table 3. 


TABLE 3 
Development time Parental ages 

Crosses (in days at 17°C) (in days) 

2e dd Le dd 
Ye/vg 23 23 | 5 5 
Y2/Og 23 26 De ae 21 
O2/Og 25 25 24 24 
OS d 26 23 22 5 
FS/F3 6 5 | 18 18 
F2/Sd 8 16 19 | 35 
S2/Sd - 4 26 26 
SP/Fd 16 8 35 20 


The theoretically ideal conditions could not be exactly realized 
for technical reasons. Results for viability and development time are 
given in table 4 and fig. 1 and 2. 

The viabilities of the same strain but with different parental ages 
and parental development times differ considerably. Moreover the 
number of individuals hatched from 300 eggs varies more in the Fi 
of crosses with different parental ages (crosses Y@/Y4, YP/OG, 
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17 25 30C* 17 25 30C* 


Fig. 1. Viability. Total number of individuals hatched from 300 eggs. 
: young; O: old; F: short development time; S: long development time. 
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17 25 30 C° 17 25 30C* 


Fig. 2. Development time of the females. Symbolism as in Diem. 
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TABLE 4 
ne es 
Temperature Temperature 

Crosses Crosses 

17e 25° 30°C 1e 20 30°C 
CON 214 175 ame ROES 216 170 119 
Yo/Og 233 240 115 F2/S& 214 Zil 89 
OL/Og 274 212 204 SPISd 217 204 147 
ONT 149 146 77 SQIF4 195 174 151 
Yo/Ydg 2202 11.6 10.8 FO/F4 28.5 12.8 10.9 
Ye/O8 23:5 12.6 WE FOIS 23.5 13.4 10.5 
OL/Og Zl Lis 10.0 SP/Sd 32.5 LSA 22 
OL/Y 4 AAE hl 8.8 SQIF4 21.6 10.0 9.0 


O2P/Og, OP/YJ) than in the F; of crosses where the parental ages 
are identical and where only the parental development time varies 
(crosses FOQ/F 3, FO/S4, SLQ/SH, SL/IF A). 

Filial development time, like viability, for a given strain but with 
different parental ages and parental development time, differs with 
respect to these values, with respect to temperature of course, and 
with respect to the nature of the cross within a single group (e.g. 
F9/S4 with respect to FQ/F4). In the Fi of crosses of the equal 
parental development time groups (crosses YP/Y 3, YL/O4, OLP/OF, 
OP/Y 4) the differences between the development times are negligible 
in relation to the differences in the unequal parental development 
times group (crosses FQ/F 3, FQ/S4, SLQ/S3, SL/F 4) which range from 
about 3 days at a temperature of 30°C to approximately 11 days at 
a temperature of 17°C. 

As such our experiments show an influence of the two variables, 
parental age and development time on the viabilities and the develop- 
ment times of the offspring. Consequently, in all the experiments to 
be described both variables were made as homogeneous as possible 
(see methods). Their influence, though manifest, but having been 
made constant, will be left out of consideration in the following 
pages. It should however be noted that our results show some indi- 
cations as to the existence of “balance effects’ in the developmental 
time from one generation to another as DURRANT (1955) has found 
for fecundity in Drosophila (see the SQ/FG and F@/S4 crosses, fig. 2). 
On the other hand the unadequate number of our experiments do not 
offer indication as to progressive increase or decrease, or cyclic 
5 


Genetica 
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changes of the viabilities of the offspring with relation to the parental 


ages. 


2. Temperature Conditions 


From table 5 it appears clearly that viability shows differences in 
behaviour with respect to genotype and temperature. 

_ A Canton-Inbred cross was reared with the wild strain Canton-S 
and the 25th brother-sister generation of the Inbred strain. The 
hybrid I®/C4' is constantly superior to the two homozygotes C9/C5 
and I®/I4. 


TABLE 5. Variation of viability as number of individuals hatched from 300 eggs; 
size as function of temperature in different families (explanation in text). 


Viability Size (in mms) 


Ll 25 30° 


ize | 25° | 30° 


Crosses | 
Canton-Inbved. | 


[ 
CPC 174 204 | 67 1.030 0.985 0.912 
C2/Id 150 99 128 1.010 0.920 0.936 
I9/Id 197 216 70 0.955 0.930 0.922 
I9/C4 201 254 112 1.018 0.969 0.912 1) 


Watou-Swedish. | 
We/Wg 112 116 4 0.903 | 0.890 — 


Wels 143 167 13 0.902 0.878 0.875 
SQ/S4 202 196 100 0.875 0.875 0.830 
SQ/Wd 237 92 49 0.821 0.881 0.835 2) 
Inbrved-e; br. | 

I9/Ig 90 jier 114 0.938 0.852 0.845 
I®/e; brg AIS me 88 75 0.906 0.893 0.863 

e; br®@/e; brd 149 183 o3 0.855 0.820 0.759 

e; br®/I4 60 111 41 0.961 0.887 0.889 2) 3) 
Inbved-rìi, e. 

Ie/Ig _ 90 127 114 0.938 0.852 0.845 
I®/ri, ed 111 128 154 0.889 0.887 0.756 

ri, eQ/ri, ed 164 131 31 0.899 0.870 0.850 

ri, e@/I4 176 his 43 0.887 0.818 0.794 2) 3) 


1) Size of the males. 
2) Size of the females. 
3) 17°C is replaced by 19°C. 
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— A Watou-Swedish cross was reared with the two wild strains Watou 
and Swedish-B. The hybrid W@/S4 is constantly intermediate in 
relation to the two homozygotes. 

— An Inbred-e; br cross was reared with the ebony; brown strain and 
the 30th brother-sister generation of the Inbred stock. The hybrid 
e;br@/Ig is constantly inferior to the homozygotes e;br@/e; br4 
and 19/14. 

These relatively simple relations are not common and in general the 
behaviour of the hybrids in relation to the homozygotes is much 
more complex and, at first sight, does not present any comprehensible 
reaction norm in relation to the homozygotes. For instance: 

— In the Watou-Swedish cross, the hybrid S9/W 4 is, in relation to the 
homozygotes, intermediate at 30°C, inferior at 25°C and superior at 
176: 

— In the Inbred-ri,e cross, the hybrid I®/ri,eg is, in relation ot the 
homozygotes, superior at 30°C, inferior at 25°C, intermediate at 
BC 
Like the viability results, the size results show differences in be- 

haviour. 

— In the Inbred-ri,e cross, the hybrid ri,e®@/I4 is constantly inferior 
to the ri, e@/ri,eg and 19/14 homozygotes. 

— In the Watou-Swedish cross, the hybrid W@/S4 is very probably 
constantly intermediate between the W@/W4 and S2/S4 homo- 
zygotes. 

— In the Inbred-e;br cross the hybrid e; br®/I4 is constantly superior 
to the homozygotes. 

Complex results are more common, as for instance: 

— In the Canton-Inbred cross, the hybrid is, with respect to the homo- 
zygotes, superior at 30°C, inferior at 25°C, intermediate at 17°C. 
From these results we can conclude that 1°. different genotypes, 

homogenized as to parental ages and parental development times, 

react differently, as far as their size and their viability is concerned, 
throughout a given temperature range; 2°. a given genotype has 
erratic reaction norms as a function of the variations of the environmental 
conditions (temperature) as far as its size and viabihty are concerned. 

This statement is valid for pure strains as well as for hybrids. The 

behaviour of the hybrids is in no apparent way distinct from the 

behaviour of the pure strains. 
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3. Reciprocal crosses 


Generally speaking reciprocal crosses show no differences within 
the optimal temperature range. 

— In the Inbred-e;br cross the two hybrids, at 25°C — which is the 
optimal development temperature for Drosophila melanogaster — 
for size and compared with the homozygotes ones, show an equiva- 
lent superiority (table 5). 

— In the Watou-Swedish cross the two hybrids, at 25°C too, are 
intermediate between the two homozygotes values, and are practi- 
cally identical. 

— In the Inbred-ri,e cross the viability of the two hybrids at 25°C 
does not differ significantly. 

On another side, in extreme temperature conditions, it is seen that 
reciprocal hybrids always differ, though they may be equivalent at 25°C. 
— In the Inbred-e; br cross the sizes of the reciprocal hvbrids I®Q/e; br3, 

e;br®@/I4 are almost equal at 25°C but exactly opposed at 19° C and 

SOC 
— In the Watou-Swedish cross there are the same equivalence of the 

male sizes of the reciprocal hybrids W2@Q/S4 and S9/W4 at 25°C, 

and the same opposition at 17° and 30°C. 

— In the Inbred-ri,e cross we found the same behaviour in the 
viabilities: equivalence at 25°C, complete opposition at the extreme 
temperatures. 

Finally, in other crosses the reciprocal hybrids do show differences 
even under optimal temperature conditions. 

— In the Canton-Inbred cross the reciprocal hybrids C4/14, I9/C4 are 
almost identical at the extreme temperatures, they are opposed at 
the 25°C temperature, i.e. for the females size. 

— The same behaviour occurs in the Inbred-ri,e cross for the size of 
the males. Approximate equivalence of the hybrids at the extreme 
temperatures, opposition at 25°C. 

— The viability of the reciprocal hybrids in the Swedish-Watou cross 
is identical at the superior extreme temperature, opposed at 25°C, 
and at 17°C the positions are reversed. 

Summarizing briefly, it can be said that reciprocal hybrids always 
differ in their reaction norm with respect to every character studied, i.e. 
size and viability, and in the given vange of temperatures. 
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4. Preliminary Discussion and Conclusions 


Momentarily our purpose is not to analyse the very intricate 
reactions of Drosophila melanogaster homozygotes and reciprocal 
hybrids on the large scale of. temperatures used, nor to explain the, 
at first sight, unintelligible and unruly responses of a given nucleo- 
cytoplasmic combination to systematic changes of external environ- 
mental conditions, nor to account for the different reactions of a 
given genome in different cytoplasms. That will be done later. For 
the moment the important features are following : 

1®. The viability and the developpent time of a zygote depend 
respectively and at least partially on the parental ages and parental 
development times. Therefore and in order to avoid observational 
bias, parental ages and parental development times of the strains 
used should be taken into account. This is to be sure the less im- 
portant fact because they may be taken into account unconsciously 
for practical purposes and technical readiness. It must be noted that 
DuRRANT (1955) concludes to similar facts in an opposite manner, 
asking for better sampling technique. Since a systematic study of 
the effets of parental ages and development times was never under- 
taken, and since the reaction norm of the offspring seems not to be a 
constant as a function of parental ages and development times, that 
conclusion, introducing a new unknown variable in genetic studies, 
seems to be unjustified. 

2°. Although the incredible lack of information concerning differ- 
ences between reciprocal crosses, some reciprocal hybrid differences 
have been found in mice (MARSHAK 1956, Mac LAREN and MicHIe 
1954, 1955, 1956), in Drosopmla melanogaster (Nrcorerri and 
SorIMA 1956), in Drosopmila subobscura (CLARK and SMITH 1955), 
in Epilobiwm (MicHarris 1943, 1951, 1954) in Galeopsis (HAGBERG 
1952a and b). That lack of information is easily understandable if we 
underline the fact that, on the one side, in many genetic experiments 
the reciprocal crosses are neglected, and that, on the other side, in 
many genetic experiments the optimal environment is the only one 
used. As a matter of fact, it is known already that extreme conditions 
favour the exteriorization of differences between reciprocal hybrids, 
for instance, high temperatures, long days and intense illumination 
(FURTAVER 1940), insufficient manuring and low temperatures com- 
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bined with short days (MicHarrIs 1935, 1939). Extreme climatic 
conditions have the same effect in interspecific Lycopersicon hybrids 
(Lewis 1952). This fact is confirmed for intraspecific reciprocal hybrids 
in our experiments, considering the large differences of reciprocal hy- 
brids in both extreme temperature combined with severe competition. 
In the medium temperature range this last factor acts alone and may 
be responsible for reciprocal differences at the optimal temperature. 

It seems justified to conclude that the cytoplasmic origin de- 
termines the differences of reaction norms in reciprocal hybrids, 
which generally are not made or not separated by many authors, but 
that such differences of reaction norms are only observable in ap- 
propriate, extreme external conditions. We can thus agree with 
Heuts’ (1956) conclusions based on his study on the reactions of 
racial hybrids of Gasterosteus aculeatus in different salinity — temper- 
ature conditions: “Considering the fact that experimentation in a 
single set of environmental conditions is rather usual practise in 
genetics, our data throw justifiable doubts on the validity of most, if 
not all, claims concerning the absence of cytoplasmic differences in 
crossing experiments’ (Heuts 1956). 

3°. A non-constant behaviour of hybrids with respect to homozygous 
parents have been found several times: in Drosopmila subobscura in 
relation to temperature (CLARK and SMITH 1955), in barley in relation 
to density of culture and temperature (GUSTAFSSON, NyYBOM and Von 
WETTSTEIN 1950; NyBoM 1950, GUSTAFSSON and NyBoM 1950). It is 
quite clear that we find here the same situation as already quoted in 
our Drosopmla experiments, that is to say, a particular nucleo- 
cytoplasmic combination, homozygote or not, which, following the 
environmental conditions, changes its reaction norm in relation to 
another, homozygote or not. LERNER (1954) has summarized in his 
book “Genetic Homeostasis”’ the evidence concerning such phenomena, 
and has concluded to a general higher independance of heterozygotes 
on external conditions, compared to homozygotes. Our results have 
an obvious bearing on LERNER’s hypothesis. This bearing will be 
discussed later. Suffice it to state that our data do not corroborate 
directly this hypothesis. 

4°. It is clear, a priori, that no simple and verbal explanation is 
applicable to the facts disclosed. As an example of mere verbal 
explanation and tautological description of the phenomena like those 
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observed by us may be taken the terms “overdominance”’ or “‘super- 
dominance”’ and “superrecessivity”’. This notions were introduced 
by Hurr (1945) to cover the heterosis phenomenon, supposedly due 
only to allelic interactions. The lack of precision of such terms, as 
well as their incapacity to describe gene-cytoplasm interactions 
which our results disclose is a sufficient reason to discard them. 
A hypothesis covering the evidence for which that “superdominance”’ 
notion was applied, as well as the nucleo-cytoplasmic interactions 
will be discussed below. 


IT. Explanatory hypothesis of the phenomenon. 
Nucleo-cytoplasmic and genetic interactions. 


To the very intricate situation, described in the previous paragraph 
and concerning the variations in behaviour of an hybrid with regard 
to environmental conditions and cytoplasmic origin, we have applied 
an hypothesis of Heuts (see SPAAS and Heuts 1958), by which 
viability can be derived, in Salmonidae and their hybrids, from a 
differential interaction between certain nuclear and cytoplasmic 
properties at different temperatures. 


L. Viabilsty 
A. The +(+), eele), + +(e) and eel +) strains 


(+) represents an individual of ++ genomic formula, Canton- 
S, with a (+) cytoplasm, Canton S. The same applies for ee(e). On 
the other hand in the +-—-(e) strain half the J-+ genome was 
introduced in a foreign (e) cytoplasm, and vice-versa in the ee(+) 
strain half the ee genome in the foreign (+) cytoplasm. 

The development speeds of the two homozygotes +(+) and 
ee(e) strain are given in table 6. 


TABLE 6. Development speeds X roo of the +(+) 
and eele) strains. Males and females 


Strains dd oe 
ER HNO 7 3.49 3.53 
25 6.67 6.75 
30° 8.24 8.54 
TTT 4.14 4.34 
20 7.99 8.11 


30° IE) 10.14 
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The viability of the four strains (fig. 4a for the females, 6a for the 
males) are given in table 7. 
TABLE 7. Viability — number of hatched flies, males 


and females, from 300 eggs — of the +(+), eele), 
dt (e) and eel) strains 


Strains dd 2e 
atten 18 101 110 
25u 104 118 

SO 86 86 

ee(e) (Zie 46 sl 
Zon 38 41 

30° 55) dd 
selen 108 85 
Zok 92 130 

30° 75 105 

eel) ijt 71 104 
255 75 88 

30° 56 74 


The logarithms of the speeds of development of the females of the 
two pure nucleo-cytoplasmic combinations +—+(—) and ee(e) (table 6 
— fig. 3) are plotted against the inverse of the absolute temperatures. 

Thereafter, considering only the real development speeds of the 
female +(+), and by visual adjustement, a straight line (continuous 
line on fig. 3) is traced, which deviates, at the ordinate values 32.98 — 


PIN x 10.000), 33.53 (25°) and 34.45 (17°), from the real 
30° + 273.2 


development speed values of the females. When the differences be- 
tween the values of the straight line and the real development speed 
values are expressed as a percentage of the values on that line, that 
percentage subtracted from 100 and transferred on a graph, a curve 
is obtained which is represented by a thick continuous line (fig. 4b). 
It is possible to adjust a second straight line (interrupted line on 
fig. 3) to the development rates of ee(e). The differences between the 
values of that second straight line and the real development speed 
values of ee(e) (closed circles) expressed as a percentage of the values 
of that straight line, that percentage subtracted from 100 and 


recorded on the same graph (fig. 4b) gives a curve represented by a 
fine continuous line, 
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If we return to the first straight line (continuous line on fig. 3) 
and ignore the differences between the values of that straight line and 
the real development speed values of + +(+) but consider the differ- 
ences between that straight line and the values of ee(e), which we 
expressed as a percentage of the values of the straight line, and if we 
substract that percentage from 100, we obtain (fig. 4b) a third curve 
represented by a thick interrupted line. 


speed 


log. development 
Ed 


32,98 3353 3445 1/T° 


Fig. 3. Continuous line: ++ “theoretical'’ development speeds. Interrupted 
line: ee “theoretical development speeds. Open circles: real development 
speeds of +(+). Closed circles: real development speeds of ee(e). Females. 
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Fig. 4a and b. Observed and calculated viabilities of the +(+), +—+(e), 
ee(e) and ee(+) strains. Females. 
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Finally from the second adjusted line, the differences between the 
latter and the values of + +(+) expressed as before, give a last curve 
(fig. 4b) represented by a fine interrupted line. 


TABLE 8 


— Real development speed values of the strains ++ (+), ee(e). 

— Values of the two calculated straight lines. (Values of the theoretical development 
speeds). 

— Differences expressed im percentage of the values of the calculated straight line 


(mortality). 
— Differences in percentage subtracted of Too (viability). 
— Females. 
Real Calculated Theoretical Theoretical 
Strain development | development mortality viability 
speeds speeds 98 96 
H(A) 17° 3.53 3.80 Zal O2 
25° 6.75 7.50 10.0 90.0 
308 8.54 11.30 24.2 75.8 
ee(e) lan 4,34 3.25 6160) 66.5 
200 8.11 5.15 57.4 42.6 
30° 10.15 6.80 49.1 50.9 
Jt(e) 17° — — 14.2 85.8 
25e — — 8.1 ee 
30° En — 10.2 89.8 
eelt) la — — 8.6 91,4 
202 — — 31.0 69.0 
30° — — 25:5 74.5 


Graph 4b is a verv similar to graph 4a. 

This means that the differences between the first adjusted straight 
line (continuous line) and the values, at each temperature, of the real 
development speeds of + +(+) represent the mortality of the + +(+) 
strain. Likewise the differences between the second straight line (inter- 
rupted line) and the real development speed values of ee(e) represent 
the mortality, at each temperature, of the ee(e) strain. 

Morcover the differences between the first straight line and the real 
development speed values of ee(e) represent the mortality of the 
+ +(e) strain; the differences between the second and + + (+) repre- 
sent the mortality of the ee(e) strain. It appears thus that the different 
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viabilities of two reciprocal nucleo-cytoplasmic hybrids [4 H(+) and 

+ (e) on the one hand; ee(e) and ee(+) on the other] can each be 

derived from the deviations of the developmental rates of their 

respective nucleo-cytoplasmic combination of which the origin of 

cytoplasm and nucleus is identical (+ +(+) and ee(e) on the one hand: 

ee(e) and +(+) on the other] from a single straight function of the 

absolute temperatures (fig. 3, continuous line on the one hand; inter- 

rupted line on the other). Let us note that there is one and only one 

possible straight line the differences of which, at each temperature, with 

the real development speeds values of ++-(+) and ee(e), are repre- 

sentative of the mortality of the +(+), ee(), and ee(e), ++ (e) 

strains. 

The following hypotheses seem to be reasonable: 

1°. The real development speeds are, at each temperature, the 
resultants of two development speeds. 

2°. A first series of these development speeds, expressed by its 
logarithmic values, constitutes a linear function of the inverse 
of the absolute temperatures, which means that these develop- 
ment speeds are submitted to a constant acceleration by the 
absolute temperature. 

3°. A second series of development speeds deviates according to the 
temperatures, on either side of the first series. The extent of 
taht deviation is unknown, but the deviation of the real develop- 
ment speeds in relation to the first series of development speeds 
is a simple function of that second series of development 
speeds. 

4°. From the fact that the viabilities of two reciprocal nucleo- 
cytoplasmic hybrids, differing only by the cytoplasm, can be 
derived from a single straight line suggest that this line repre- 
sents some reaction of the developmental rates to absolute 
temperature dependent on the nuclear elements and probably on 
the chromosomes. We call it henceforth the “nuclear acceleration 
Vine’ without claiming to know the concrete mieaning of that 
expression. 

5°. The fact that the developmental rates of two reciprocal nucleo- 
cytoplasmic hybrids, differing only by the cytoplasm, are never- 
theless different, and that their equally different viabilities can 
be derived from a single nuclear acceleration line and from the 
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two different series of maternal development speeds, suggest 
strongly that the observed viabilities are also influenced by the 
cytoplasm. A reasonable hypothesis would seem to be that the 
second series of speeds is only transmitted by the cytoplasm. 
That second series will be called “cytoplasmic acceleration curve” 
again without claiming to know the concrete meaning of that 
expression. Thus, as a result of this first experiment we should 
admit that the ee(e) strain differs from the +(+) strain not 
only by the e gene, but also by the (e) cytoplasm. 

6°. At each temperature, the differences, i.e. the greater or smaller 
distance, between the development speed of a given genome and the 
development speed of the cytoplasm into which it is introduced 
determines the higher or lower mortality of a given genetic, te. 
nucleo-cytoplasmaic combination. 


Following the same procedure and starting from the real develop- 
ment speeds of +(+) and ee(e) strains (table 9) nuclear linear 
acceleration lines (fig. 5) can also be adjusted for the males of the 
HH (H), eele), ++H(e), ee(+) strains (fig. 6a and b). The results are 
again concordant. 

TABLE 9 


— Real development speed values of the strains +(+), ee(e). 
— Values of the theoretical development speeds of the ++ and ee genomes. 
— Theoretical mortality and viability. 


— Males. 
Real Theoretical Theoretical Theoretical 
Strains development | development mortality viability 
| speeds speeds Se DA 
dE fe A 3.49 4.00 | 1207 87.3 
250 6.67 6.47 3.0 97.0 
30° 8.24 8.70 9.2 94,8 
ee(e) les 4.14 2.80 47.8 E22 
258 7.99 Ball 59.8 40.2 
30° CITO 7.40 35.0 65.0 
d-t(e) 17° — — 3.5 96.5 
20 — — 23.4 76.6 
30° — — 16.5 83.5 
eelt) IA — — 24.6 75.4 
25° — — 30.0 70.0 
30° — — 11.3 sle 


eN 


INTERACTIONS IN DROSOPHILA MELANOGASTER 207 


log. development speed 


32,98 33,53 3445 1/T° 


Fig. 5. Continuous line: ++ “theoretical” development speeds. Interrupted 
line: ee “theoretical development speeds. Open circles: real development 
speeds of +(+). Closed circles: real development speeds of ee(e). Males. 


ee 


30 


Observed Calculated 


17 25 30C* 17 25 30C° 


Fig. 6a and b. Observed and calculated viability of the +(+), ++(e), 
ee(e) and ee(+) strains. Males. 


B. The Watou-Swedish cross 


In the experiment described above the +(+) and ee(e) strains 
were different only for the e gene and, as it appeared later, for the 
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(e) cytoplasm. The Watou-Swedish cross was made from the two wild 
stocks Watou and Swedish B which, on the contrary, are supposed to 
differ by a great number of genetic constituents. 

As in the former instance we can adjust the nuclear acceleration 
line from the real development speeds of the WW(W) and SS(S) 
genetic combinations (table 10, Ist row). This nuclear acceleration 
line gives the theoretical development speeds of the WW and SS 
genome (fig. 7). The differences between the nuclear acceleration line 
and the real development speeds — a simple function of the cyto- 
plasmic acceleration curve — is truly representative of the mortality 
of the WW(W) and SS(S) strains. 

Again, a single nuclear acceleration line WS is shared by the two 


log. development speed 
EN 


32,98 33,53 34,45 

UT 
Fig. 7. “Theoretical” development speeds of the WW, SS, WS genomes. Open 
circles: real development speeds of SS(S). Closed circles: real development 
speeds of WW(W). Females. 


RE 
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reciprocal hybrids W®/S4 and S9/Wg4. On the one hand, the differ- 
ences between that straight line WS and the real development speeds 
of WW(W) are representative of the mortality of the W@/S& hybrid, 
with cytoplasm (W); on the other hand, the differences between the 
nuclear acceleration line WS and the real development speeds of 
SS(S) are suitably representative of the mortality of the SQ/W& 
hybrid with cytoplasm (S) (table 10, fig. 8a and b). 

The new fact discovered in this experiment is that the nuclear 
acceleration line of a hybrid nucleus differs considerably from both 
parental nuclear rates attributed to the homozygotes nuclei. The 
characteristics of the hybrid nuclear rate can apparently not be 
predicted by additive arithmetic from those of the homozygotes. 
The viability of a given hybrid combination seems to be determined 
by the distances between the maternal cytoplasmic acceleration curve 
and a nuclear acceleration line exhibiting new characteristics through 
interaction of the parental halve nuclei. The real maternal cytoplasmic 
acceleration curve is, of course, unknown. The real development speed 
is probably to be considered as resulting, at each temperature, from 
another interaction between the real “cytoplasmic’ developmental 
speed, and the empirically known “nuclear” rate. Inevitably, however, 


Calculated 


5 TE 3oct 17 25 30C* 


Fig. 8a and b. Observed and calculated viabilities of the WW(W), WS(W), 
SS(S) and SW(S) strains. Females. 
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TABLE 10. As in the tables 8 and 9. 
Watou-Swedish cross. Females. 


NN ( 


Real |_ Theoretical Theoretical Theoretical 
Strains development | development mortalities viabilities 
speeds speeds Do 9 
WW(W) 17° 5.40 5.00 8.0 92.0 
25 7.24 6.60 9.6 90.4 
30° 9.09 7.85 7.0 83.0 
SS(S) log 4.79 4.60 4.1 95.9 
25 6.15 6.47 4.9 95d 
30° | 9.01 18 Wee 80.9 
WS(W) 172 — 4.92 SEL 90.3 
25 — 6.85 5.6 94.4 
30° — 8.44 | 1D 92.5 
SMS) MELS — 4.92 2.6 97.4 
25° — 6.85 10.2 89.8 
30° — 8.44 12.5 87.5 


the use of the observed developmental rates instead of the unknown 
cytoplasmic rates leads to the consequence that the observed viabili- 
ties are covered by the calculated viabilities in absolute terms. 
However, the value of the hypothesis used is born out by the possi- 
bility it offers to predict the general behaviour of a given genetic 
combination as well as the precise comparative behaviour of hybrids 
with respect to the parental strains. 


C. Canton and Swedish strains. SS(C)12 and CC(S)12 crosses 


The SS(C)12 and CC(S)12 strains have been obtained from the 
repeated reintroduction during 12 generations of the half-nucleus S 
into the (C) cytoplasm, and of the half-nucleus C into the (S) cyto- 
plasm. We can thus venture that the CC nucleus is entirely replaced 
by the SS nucleus, and vice-versa the SS nucleus by the CC nucleus; 
hence the SS(C)i2 and CC(S)i2 symbolism (see also Material and 
Methods, table 2). 

The nuclear acceleration lines of the CC and SS genomes are 
adjusted as ín the two former cases (table 11, fig. 9, 10a and b). 

Assuming that the SS genome is entirely replaced by the CC genome 
and vice-versa, and assuming further that the original maternal 
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ot 


log. development speed 


32,98 3353 3445 1/T° 


Fig. 9. Theoretical development speeds of the CC, SS, “CC(S)12”’ and “SS(C)12”’ 
genomes. Open circles: real development speeds of CC(C). Closed circles: real 
development speeds of SS(S). Males. 


cytoplasmic properties have been carried autonomously through the 
successive generations the viabilities of the two series of descendents 
should be obtained by calculation based on a nuclear CC and SS 
acceleration line, and on a (S) and (C) cytoplasmic curve respectively, 
the latter being derived from the SS(S) and CC(C) developmental 
rates. Neither one of both assumptions proved to be true. We were 


no 
90 
70 


50 


ES CC(S2 
SS(S) 
ei US HCHI 


30 


Calculated 


10 Observed 


17 25 30C* 17 25 30C* 


Fig. 10a and b. Observed and calculated viabilities of the CC(C), “CC(S)12”’, 
SS(S), “SS(C)12’’ strains. Males. 
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TABLE 11. As in the tables 8, g and To. 
CC, SS, CC(S)12 and SS(C)12 crosses. Males. 


Real Theoretical Theoretical Theoretical 

Strains development | development mortalities viabilities 
speeds speeds 95 08 
SE) AE 3.69 2.80 29.6 70.4 
25 7.94 6.80 16.7 83.3 
30° 9.24 11.60 20.3 79.7 
SS(S MZ 3.23 4.67 30.8 | 69.2 
258 7.61 5.68 33.9 66. 1 
30° 8.00 6.48 23.4 76.6 
CC(S)1a 17° == | 4.10 10.0 90.0 
Zok —= 135 8.0 92.0 
30° — 10.50 12.0 88.0 
SS(C)1i2 17° == 5.15 372 62.8 
256 DT 5.68 33.9 66.1 
30° — 6.05 32.2 67.8 


unable indeed, to deduce the viabilities of the respective crosses, 
either by combination of the SS and CC nuclear acceleration with the 
developmental rates of CC(C) and SS(S), or by the combination of the 
same nuclear accelerations with respectively the SS(S) and CC(C) 
developmental rates. 

The only reliably representation of the so-called CC(S)i2 and 
SS(C)12 strains viabilities are obtained by calculation from new 
nuclear acceleration lines, specific for both crosses, and the develop- 
mental rates of CC(C) and SS(S) respectively. 

These results suggest that the repeated reintrodunction of the 
half-nucleus Canton in the Swedish cytoplasm, and conversely of the 
half-nucleus Swedish in the Canton cytoplasm have transformed the 
Swedish and Canton cytoplasms in the direction of the Canton and 
Swedish cytoplasms. It would appear therefore that the cytoplasm is 
not an autonomous entity and on the contrary that there exists for 
each nucleus one cytoplasm adapted to that nucleus. (The exact 
meaning of that “adaptation” will be discussed later). These results 
confirm the hypothesis expressed with respect to the differentiated 
cytoplasm of CC(C) in relation to +(+), ee(e) being only a mutation 
of HH (+-) however. On the other hand they show that these cyto- 
plasmic differences are controlled by the nucleus and that the cyto- 
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plasm should rather be considered as an external and labile factor in 
relation to a given genome. 

On the basis of our data the following conclusions seem to be 
justifiable: the cytoplasmic differentiation, in relation with one or 
several gene differences is a fact. These cytoplasmic differences, however, 
are not autonomous. They appear or disappear under the influence of the 
nucleus. Therefore, since the cytoplasm, in our material, lacks autono- 
mous constancy it cannot be properly called a hereditary material. 


D. Note on the influence of crowding on the expression of viability 


We have noted that the agreement between the observed and 
calculated viabilities exists only in the general behaviours of the two 
compared curves, observed and calculated, and we have given the 
reasons for it. Moreover, if the calculated ratios between the viabilities 
of the diverse nucleo-cytoplasmic combinations of a given cross are 
sufficiently representative of the observed ratios of the same combi- 
nations, there are still some calculated points — see for instance the 
calculated viability curve of the ee(+) males at 30°C (fig. 4a and b) — 
which do not agree with the reality. 

We may recall that all our experiments have been made in a highly 
crowded environment (300 eggs for a standardized surface of 5,47 cm?). 
Vital competition was certainly a factor of mortality. In fact these 
experiments should have been carried out a second time in a normal 
environment, without crowding and possible competition. 

As on one hand, we have shown that the viability of a given genetic 
combination depends on a nucleo-cytoplasmic interaction; and on the 
other hand we have shown that the cytoplasm is a labile factor 
reacting in a irregular way to regular variations of a given environ- 
ment, it seems reasonable to conclude that the nucleo-cytoplasmic 
interaction which is expressed with respect to temperature as the 
calculated viability — and, more particularly, that the cytoplasmic 
action is influenced by the highly crowded environment in which the 
experiments have been conducted. This new variable, as yet indefi- 
nite, is certainly not a constant and there probably exists for each 
observed genetic combination, and with respect to temperature, a 
positive or negative coefficient which would account for its variable 


influence. 
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We actually ignore the value of that coefficient. Its existence, 
however, can be demonstrated indirectly. Indeed, in experiments 
carried out on Salmo trutta fario eggs (SPAAS and Heuts 1958) at 
five different temperatures, where, of course, all crowding effects 
are excluded, the agreement between calculated and observed 
viabilities is perfect (table 12, fig. 11). This shows the considerable and 


TABLE 12. Real and theovetical nuclear development speeds in 
Salmo trutta fario. Theovetical movtalities. 
(after SPAAS and Heuts 1958) 


Real Theoretical Theoretical 
oC development | development mortalities 
speeds speeds 0E 
4 0.93 Öa4 6.7 
7.8 1.80 1.60 1e 
or 2.12 2.04 39 
10.6 2.56 2.60 1.5 
12.5 3.33 3:52 6.7 
% 
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End 
- ee 
Oe an Els SS 
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10 Se 2 
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fig. 11. Observed viabilities and calculated mortalities in Salmo tutta fario 
(after Spaas and Heuts 1958). 
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till now not formulable influence of crowding in our experiments !). 


2. Relation between size, fecundity and viability 


Viability seems to be precisely determined by the differences be- 
tween the nuclear acceleration line and the cytoplasmic acceleration 
curve. It is therefore fitting to see if size and fecundity are determined 
by the same factor, or in any case, how viability, size and fecundity 
are related to each other. 


A. Size — Viability 


Size is expressed by the length of the thorax, from the anterior 
end of the mesonotum to the posterior end of the scutellum. The mean 
size of the descendents of a given cross at a given temperature is 
based on the measurements of thirty individuals of each sex. For a 
given cross samples of less than ten individuals are discarded. Six to 
seven thousand individuals have been measured. 

Correlations between size and viability have been established at 
each temperature, for males and females, starting from the average 
sizes of the different strains 2). Correlation between size and viability 


TABLE 13. Correlation between size and viability 


DE | Females | Males 
17° Öl | EEn 
25° —0.12 —0.57-* 
30° — 0.02 | —0.40-- 


ARP O1 
ww OOTAP. n0:02 


1) Other results exist which are in good agreement with those given above and 
in the first place for the males of the Watou-Swedish cross and for the females 
of the CC(S)12 and SS(C)12 crosses; in the second place for the females and the 
males of Inbred-ri; e crosses and Inbred-e, br crosses. 

It should be noted that the nuclear acceleration lines for males and females 
of a given cross are not identical. 

2) Following groups of crosses were included in these and further calculations: 
Swedish-Watou, Inbred-e;br, Inbred-ri;e, Inbred-Canton, the ++ (+4)/ee(e) 
group (see under IV: comparative study of cytoplasms), the CS12/SC12 group 
(see under IV), the Young-Old series and the Slow-Fast series of crosses. 
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is negative for males and females at the three temperatures of 17, 25 
and 30°C. It is highly significant (P < 0.01) for females at 17°C, for 
males at 25° and 30°C; significant (P < 0.02) for males at 17°C; 
negative but not significant in two cases: females at 25° and 30°C. 

Viability, if one neglects all environmental conditions, such as 
density of population for instance, is basically determined by the 
nucleo-cytoplasmic formula of the cross, and more particularly by the 
differences between the nuclear acceleration line and the cytoplasmic 
acceleration curve as has been demonstrated. Viability of a given 
strain is therefore determined by the zygotic constitution of that 
strain. In the experimental conditions used with Drosophila the 
greater or lesser viability of a cross introduces a new environmental 
factor namely a greater or lesser density of population, which in turn 
leads to less or more severe competition influencing directly the sizes 
of the flies. This explains sufficiently the negative statistical relation 
between size and viability. Henceforth it seems logical, unless further 
information is forthcoming, to consider that in our experiments size 
is variable, directly dependent on competition, and indirectly on 
viability. 


B. Fecundity — Size 


Fecundity, or average daily egg-laying, was calculated according to 
two different methods, whose results are comparable and have the 
same meaning (see Materials and Methods). About 150.000 eggs have 
been counted in these experiments. 

Correlation between fecundity and size was calculated for females 
reared at 25°C — which is the optimal temperature for Drosophila 
melanogaster development — and kept for the measurement of fecundity 
at that temperature of 25°C; the correlation between fecundity and 
size is positive and highly significant (r — 0.50; P < 0.01). 

Other correlations for individuals reared at 17°, 25° and 30°C and 
observed for fecundity at temperatures other than those at which the 
development occurs have also been calculated; they are not significant, 
except in two cases. For females reared at 30°C and transferred for 
fecundity at 25°C there is a slight negative correlation (r = —0.39; 
P < 0.10). That point will be explained below (in paragraph C). 

So, at least at the optimum temperature of 25°C, size and fecundity 
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are closely related. Another fact which leads to the same conclusion 
is that the correlation between size and fecundity for females reared at 
25°C and observed for viability at 17°Cis non-significant (rtecundity.size 
— 0.16) but for a fixed viability the partial correlation is highly 
significant (rrecunaity-size.vrabiuty — 0.64; P < 0.01). Without fuller 
information, it is, however, impossible to know the dependency 
relations of the two variables, size and fecundity. 


C. Relations between size, fecundity and viability 


In paragraph A the average size of the individuals was found to be 
a dependent variable of the viability of the cross involved. Both traits 
are negatively correlated and their relationship can be expressed by a 
regression line with a negative coefficient. It seemed logical to at- 
tribute this general negative relation to crowding. Now, single families 
deviate considerably in a positive or negative direction from the 
established regression line. 

The SS(C)12P/SS(C)123 cross (fig. 12, no. 24) has a viability of 42%, 
(at 25°). According to the regression line of the size in function of 
viability, the average size of that family should be of 1.029 mm; it 
is in fact of 1.068 mm; the same is true for the C®/C4 (no. 20), 
CP/CC(S)114 (no. 6), and CC(S)11P/S& (no. 4) crosses which are all of 
a size superior to the average one. Other families have for a given 
viability a size inferior to the expected one. Such is the case for 
instance for the CQ/I4 family (fig. 12, no. 13) which has a viability of 
52%; according the regression line of size in function of viability the 
size of CQ/I family is expected to be 1.026 mm; in fact it is of 
1.010 mm. The same is true for the IQ/I4 (no. 15) and SP/CC(S)11d 
(no. 5) families. 

We can now consider the negative regression coefficient as a mean 
correction factor for crowding effects, and regard the positive and 
negative deviations from the average regression as indicating a 
positive or negative heterosis for size in each one of the different 
crowding conditions. This view supposes that we take these deviations 
as resulting not from random external conditions but from genetical 
factors. SpPAAs and Heuts (1958) have shown in trout and salmon 
development, where all competition was excluded, that size as well 
as viability are both directly dependent on the nucleo-cytoplasmic 
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Fig. 12. Scatter diagram and size-viability regression line. 


interactions. These results support our hypothesis that in Drosoplmla 
too, once the crowding effect corrected by the regression coefficient, 
the remaining variations are due to the same nucleo-cytoplasmic 
interactions. So the real size, as opposed to the average theoretical size, 
of a given cross, in a high-crowding environment, is determined both by 
the environmental conditions and by the nucleo-cytoplasmic formula of 
the strain. 


If we are right in considering the positive or negative deviations 
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from the size on viability regression as indicative of genetical differ- 
ences, these deviations as such could show in turn correlation with 
other measured traits. We know already that fecundity as such is 
highly and positively correlated to size. From this correlation two 
regressions can be calculated (fig. 13), one of fecundity on size where 
fecundity is considered as the dependent variable, another of size on 
fecundity where size is the dependent variable. The very important 
fact is now that the deviations from the fecundity on size regression 
do not show any correlation with the higher mentioned deviations 
from the size on viability regression line, while the deviations from 
the size on fecundity regression line are highly and significantly 
correlated with the deviations from the size on viability regression 
line (r = —0.46; 0.01 < P < 0.02). 

This means that the deviations of the real size around the average 
size, in a positive or negative direction parallel the variations of 
fecundity which will deviate, in a negative or positive direction; this 
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Fig. 13. Scatter diagram and fecundity-size regression line. 


220 F. A. LINTS, NUCLEO-CYTOPLASMIC 


is on the condition that we consider size as a dependent function of 
the fecundity. 

Thus if we consider the crosses SS(C)19P/SS(C)123 (fig. 13, no. 24), 
CQ/C& (no. 20), CR/CC(S)11d' (no. 6), and CC(S)uL/Sd (no. 4) in 
relation to the regression line of size as a function of fecundity we can 
immediately see that for the given sizes — which are, as we know, 
superior to the average size — the fecundity is inferior to the average 
fecundity. Likewise and inversely the C2Q/13 (fig. 13, no. 13), I9/I4 
(no. 15) and S@/CC(S)11d (no. 5) crosses although in actual size 
inferior to the average size have a fecundity superior to the average- 
fecundity. 

It follows that size depends on the following factors: the environment, 
in this case competition, and the genetical factors, that is the nucleo- 
cytoplasmic formula, or more precisely, the differences between the 
nuclear and the cytoplasmic rates, and finally the fecundity. Thús means 
further that the traits viability, size and fecundity are all related with 
each other if, and only 1f, size is considered as a dependent variable on 
fecundity. These results mean finally that in each case considered, 
the deviations from the average regressions cannot be due to chance. 
Furthermore it would be extremely difficult to explain our data on 
the basis of environmental effects. An environmental factor inducing 
simultaneously a deviation in a positive direction for size and in a 
negative direction for fecundity, and vice-versa, would be barely 
intelligible. 

Two other results should be particularly stressed: 1°. The regression 
coefficient of size on viability is significantly higher (b — — 1.292) 
than the regression coefficient of fecundity on viability (b — 0.213); 
2°. The standard deviation of the regression line of size on viability 
is much higher (S —= 25.79) than the standard deviation of the 
regression line of fecundity on viability (S — 10.01). 

This means that: 

1°. In the experimental condition used, and for a given nucleo- 
cytoplasmic formula, it is the viability of that combination which 
determines directly the higher or lower density of population, i.e. the 
environmental conditions; so the regression coefficient of the size or 
of the fecundity on viability is representative of the higher or lower 
dependance of size and fecundity on the environmental conditions; 
the regression coefficient of size on viability being higher than 
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that of fecundity on viability it must be concluded that: in the 
high-crowding experimental conditions used, and in opposite drrections, 
size is dependent upon the environmental conditions to a greater extent 
than fecundity. 

2°. The standard deviation of the regression line of size and fe- 
cundity on viability is representative of the dispersion of the real 
values of size and fecundity. Since we have established that the 
deviations of the real size around the average size, in a positive or 
negative direction, are reflected in the expression of fecundity which 
varies around the average fecundity in a negative or positive direction, 
these deviations are not the result of chance. Thus these deviations 
can only be the result of zygotic factors, or more precisely, of the 
differences between the nuclear acceleration line and the cytoplasmic 
acceleration curve. Since the standard deviation of the regression line 
of size on viability is higher than the standard deviation of fecundity 
on viability it is clear that, in the high-crowding experimental con- 
ditions used, and in opposite directions, size is more dependent than 
fecundity on the differences between the nuclear acceleration line and 
the cvtoplasmic acceleratton curve. 

We believe that this is the correct way to explain the exceptional 
and slight negative correlation (r = —0.39; P < 0.10) between fe- 
cundity and size that exists for females grown up at 30°C and 
transferred for fecundity measurement at 25°C (see under II, 2, B). 
30°C is an upper extreme temperature for Drosophila, 31° to 32°C 
being the lethal one, and it is quite understandable that, since fe- 
cundity appears as a more constant value than size and thus as a 
more intrinsic property of an organism, there exist different thresholds, 
more rapidly attained by size than by fecundity, beyond which some 
“buffering properties’ of the zygote, namely the differences between 
the nuclear acceleration line and the cytoplasmic acceleration curve, 
are broken down so that environmental conditions are responsible 
for all or greater part of resulting size and fecundity. So, on extreme 
points of the temperature scale, the positive correlation between 
fecundity and size can be reversed. 


D. Conclusions 


The statistical treatment of our data exposed in the previous 
paragraph allow the following statements: 
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|. Viability is dependent, on each temperature, on a particular 
nucleo-cytoplasmic relationship directly determined by that temper- 
ature by differential acceleration. It is inferred that, in the experi- 
mental conditions used, the resulting viabilities, on their turn, 
constitute a varying environmental factor, by the degree of crowding 
they provoke. 

2. In opposite ways size and fecundity are grossly correlated to 
viability, i.e. crowding, the sizes diminishing in higher densities and 
vice-versa. Fecundity is less sensitive than size to variations in these 
environmental conditions. 

3. The residual variations of size and fecundity are not random. If 
size is taken as the dependent variable they are highly and negatively 
correlated with each other on the one hand; they are also correlated, 
inversely, with the zygotic factors, i.e., the nucleo-cytoplasmic re- 
lationship on the other. In the experimental conditions used fecundity 
again appears to be more independent than size. 


II. Comparative study of the genomes 


On the basis of the behaviour of reciprocal hybrids we have at- 
tributed an empirical linear temperature-rate relation to a reaction 
norm of the nucleus. A comparison of the slopes of the linear temiper- 
ature-rate relation of different genomes would seem a promising line of 
research. In view of the relative high number of different genomes used 
in this study a certain number of conclusions can already be reached. 

It is remembered that the slope of the nuclear acceleration line 
can be expressed by the Arrhenius’ constant, for which the formula is: 


log ka — log k 
red gt gp S Endais Sla, 
UT Sails 


1. Differences between sexes 


As one can see in table 14 the Arrhenius’ constants are different 
for males and females. There is, however, no regularity in these differ- 


ences. The range of variation in the constants seems to be lower for 
females. 
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TABLE 14. Constants of Arrhenius 


a. Homozygotes 


em 


2e Constant dd Constant 
WW 6.127 SS 1) 4,488 
SS 1) 7.116 | WW 6.565 
SS 2) 7.484 SS 2) 7.107 
ee 10.032 ri, e 7.268 
II 11225 e: br 8.326 
ri,e 13.662 — 10.557 
e;br 13.703 ee 13.206 
Le 14.807 II 14.273 
GE 17.429 GE 19.315 

b. Hybrids 

oe Constant dd Constant 
SS(C)12 5.202 SS(C)12 2.185 
SW 7.645 e; br-I 1553 
CG(S)12 7.866 SW 11.270 
e;br-1 10.266 CC(S)12 12.778 
ri, e-I 12.456 ri, e-I 13.234 


1) From the Canton-Swedish experiment. 
2) From the Swedish-Watou experiment. 


2. Differences between homozygotes and hvbrids 


A comparison of the homozygote strains, and more especially of the 
female values shows that the increase of temperature dependance, as 
expressed by the Arrhenius’ constants, follows roughly the sequence: 
wild stocks, mutational stocks and the CC stocks. The lowest Arrhenius 
constant is found in the Watou wild stock which was used for experi- 
ments immediately after capture, and so has not been submitted, as 
the other “wild stocks’ have been for many years, to a constant 
temperature and to a type of reproduction leading almost certainly to 
a higher degree of homozygosity than in natural populations. On the 
other hand, the upper values are those of the CC and ++ genomes, 
Le. the Canton-S stock and another derived from it by outcrossing 
to the ee stock (see Methods). It should be remembered that Canton-S 
was made isogenic (BRIDGES 1944). It seems thus that a higher 
homozygosity of the genome leads to a higher dependence of a part 
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of the zygotic constituents of temperature variations, which is 
expressed in a higher value of the constant of Arrhenius. However, 
the relatively low Arrhenius’ constant in the II strain is in every 
case lower than in the CC strain, from which it was derived by 
brother-sister mating. It seems as yet difficult to bring this fact in 
line with an hypothesis of a simple relationship between the value of 
the Arrhenius’ constant and the degree of homozygosity. 

The values of this constant obtained in hybrid strains, on the other 
hand, confirm again this view. Indeed they generally have constants 
inferior in value to those of the parental strains. That is the case for 
SS(C)12, e;br-I, ri,e-I. The genome of the CC(S)12 has about the same 
value as its parent with the lowest value. For the SW genome it is 
slightly higher than its parents; vet it should be noted that it is an 
hybrid between two wild stocks, Watou and Swedish, where Watou 
must be taken in a very special account; indeed these two wild 
stocks possess the lowest established values of the Arrhenius’ constants. 

Though it seems, according to our results, that there exists some 
relationship between the degree of homozygosity and the value of the 
Arrhenius’ constant, meaning that the rate of an unknown nuclear 
function is more independent on temperature in hybrids, it should be 
stressed that this property of the hybrids does not allow any pre- 
sumption as to their heterotic superiority with respect to viability, 
size or fecundity at the different temperatures studied. Indeed there 
is not a single significant correlation between the values of the 
Arrhenius’ constant and these properties. This is understandable if 
one bears in mind that viability, size and fecundity are not related 
directly to the “nuclear speed’’ of a genome at a given temperature, 
but to the relation between this value and the “cytoplasmic speed” 
for which the actual developmental rate is a measure. So the genome 
of the CC(S)i2 strain, males, possesses a value of the Arrhenius’ 
constant superior to that of the SS genome (table 14), and very 
inferior to that of the CC genome; nevertheless its viability (fig. 10a) 
is significantly superior, at the three temperatures of 17°, 25° and 
30°C, not only to that of CC(C) but also to that of SS(S). Therefore 
it is possible to say that higher homozygosity signifies higher temper- 
ature dependence for the genomes, but no conclusion can be reached 
about the “adaptive value” of that higher or lower dependence. 
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IV. Comparative study of cytoplasms. 


Nuclear action on the cytoplasm. 


In paragraph II, 1, A the conclusion was drawn that the cytoplasm 
of Canton-S is different from that of the ebony stock as well as from 
that of Swedish-B. 

The fate of these cytoplasmic divergences were, however, completely 
different after twelve backecross generations in both cases. Clearcut is 
the result in the case of CC(C) x SS(S) backcrosses. The technique 
used (table 2) namely repeated crossing of the CS hybrids by a male 
derived each time from one of the pure stocks has led presumably to a 
complete replacement of the genome of one of the stocks by that of 
the other, the cytoplasm being derived from the P-females. It can be 
assumed consequently that in fact CSj2 — SS and SCi2 = CC. Our 
results show that the cytoplasm as well has been transformed ac- 
cordingly, so that CS12(C) = SS(S) and SC1a(S) = CC(C). The case of 
the ee(e) and +(+) is not as clear. The backcrossing technique, 
according to table 1, results, not in a replacement of the ee-genome 
(which may be distinct for many other genes from the ++ genome) 
by the ++ genome, but in a gradual homogeneization of the existing 
differences between the original genomes in a new one. In view of the 
results obtained in the Canton-Swedish cross after 12 backcross 
generations it seems highly probable that in the case of the + +(+)/ 
ee(e) cross, and after 12 backcross generations the reciprocals only 
differ by the e gene, the remaining original differences of the ++ and 
ee genomes having been homogeneized during the backcrossing 
procedure. But since mass-mating was allowed, a certain degree of 
heterozygosity, depending on the original divergence between both 
stocks, is presumably present in each reciprocal strain. 

Since the number of backcross generations is the same in the two 
cases, and that in one case no change appears in the cytoplasm (or, to 
be exact: in the experiments carried out and with the techniques used 
no changes are perceptible), while in the other the cytoplasm is 
almost entirely transformed, it seems clear that in the one case 
[H+H(e)/ee(+)] the difference between the ++ and ee genomes 
is so small that no changes can occur in the (e) and (+) cytoplasms 
in so limited a time —it should be remembered that there is a difference 
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between (+) and (e): see II, 1, A — while in the other case (CS12/SC12) 
the differences between the CC and SS genomes are much higher and 
changes can occur in the (S) and (C) cytoplasms under the influence 


of these genomes. 
DISCUSSION 


The study of a number of important traits of Drosophila melano- 
gaster, playing a paramount role in what is generally called “fitness” 
or “selective value’, has yielded a body of results, the explanation 
of which seems impossible without an appeal to the maternal cytoplasm 
as well as to the nucleus. Besides their general implications, which 
will be discussed first, our results have obvious bearings on the 
theories of heterosis, bring new facts concerning the respective roles 
of the cytoplasm versus the nucleus, and finally permit us to gain an 
insight into the causal links between the different traits related to 
fitness. 

Appropriate reciprocal crosses and the use of an array of different 
external conditions have shown indiscutably the intervention of both 
causal complexes, maternal cytoplasm and nucleus, in the traits 
studied. Indeed viability and size in homozygotes as well as in re- 
ciprocal cross products show a seemingly arbitrary reaction, which 
cannot possibly be predicted on a classic basis of additive or epistatic 
gene effects, without recourse to a considerable number of secundary 
hypotheses. A way to elude such difficulties is the traditional one: to 
work at a single temperature and to avoid the comparison between 
reciprocal crosses. Although some authors have refused to consider 
all genetic explanation not based on reciprocal crosses in different 
external conditions (Heuts 1956, for racial hybrids) and although 
LERNER (1954, 1958) deplored the lack of genetical data for homo- 
zygotes and heterozygotes in different sets of environmental conditions, 
this tradition does not seem to be abandonned as easily. 

The systematic study of reciprocal cross products and the system- 
atic use of a wide range of temperatures against which parental and 
hybrid strain were allowed to react, resulted in a highly intricate 
picture of developmental speeds, viabilities, individual sizes and 
fecundities in the array of genotypes used. 

Now it appeared that the mortality of a given nucleo-cytoplasmic 
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combination of which cytoplasm and nucleus are of like origin }), is 
determined at each temperature by a relationship between two differ- 
ent logarithmic rates. A first appeared to be dependent on the nucleus 
and shows a constant acceleration against the inverse of absolute 
temperature. The second is dependent on the cytoplasm and describes 
an irregular curve around the “nuclear” R/T function. 

The study of first generation hybrids shows that the enormous 
differences of viabilities in reciprocal hybrids, just as in Salmonidae 
(SPAAs and Heuts 1958), can be derived, for both hybrids, from the 
different maternal rates on the one side, and from their respective 
deviations from a single nuclear function, alike for both hybrids, but 
novel with respect to those of the parents, and in each case, by no 
means predictable from the parental ones on the basis of simple 
addition. 

The nuclear rates, for pure strains as well as for hybrids, can be 
expressed by means of the constant of Arrhenius. It will be noted 
further that the exact values of the cytoplasmic rate are not known 
yet. This is essentially due to the fact that an insufficient number of 
points of the cytoplasmic curve are known by lack of sufficient 
constant temperature rooms. The values of the cytoplasmic rate are 
replaced in our calculation, for the hybrids by the real developmental 
values of the maternal parent, for the nucleo-cytoplasmic combi- 
nations with cytoplasm and nucleus of different origin by the real 
developmental values of the pure strain with the same cytoplasm as 
in the combination studied. 

We have advanced arguments to show that the developmental rate 
at a given temperature is itself the result of the interaction of two 
different and independent systems of response to the temperature 
changes, whereof the one is nuclear and the other cytoplasmic. As 
viability and development speed, individual size and probably fe- 


1) We prefer that paraphrase to the more simple notation of homozygote 
which is not without ambiguity. According to BATESON (1902) — on that matter 
see DARLINGTON and MATHER 1952 — a homozygote is a zygote derived from 
the union of gametes identical in respect of the quality, quantity and ar- 
rangement of their genes. Now a zygote of which the two chromosomal haploïd 
sets are identical, but of which the cytoplasm is of a different origin is called 
“homozygote” while ‘heterozygote’ would be more adequate. If this is 
admitted the term “hybrid” would define an organism of which one haploid 
set comes from one origin, the other and the cytoplasm being from an other one. 


Genetica 


228 F. A. LINTS, NUCLEO-CYTOPLASMIC 


cundity as well are further shown, by means of correlations, to be 
directly or indirectly related to these hypothetical system of nucleo- 
cytoplasmic interaction. Adequate experiments have revealed further- 
more that between any laboratory strain cytoplasmic differences of 
the nature described, and, of course, nuclear differences, exist. 

Our results have a first important bearing on the theories of 
heterosis, which in the recent past have known a revival of interest. 
The general characteristics of heterosis from the vast amount of 
literature on the subject can be summed up as follows: hybridization 
leads an organism to exceed both parental homozygotes in a positive 
or negative direction, on the morphological level, in size, weight, 
speed of development, fecundity, or, on the physiological level, in 
adaptability, ie. in the property of physiological homeostasis, in 
efficiency, or in the faculty to synthetize greater quantities of raw 
material from the given energetic resources (there is a good review in 
WHALEY 1944. See also GOWEN 1952). These positive or negative 
moditications have been theoretically expressed by notions called 
“genetic inertia’’ (DARLINGTON and MATHER 1949), “canalization’’ 
(WADDINGION 1953) “flexibility” (THopAy 1955), and stoutly de- 
fended by LERNER in his book “Genetic homeostasis’’ (1954). 

There is, first of all, the fact that if any homeostatic properties 
characterize the hybrids at all for as far as their viabilities, develop- 
mental speeds, individual sizes and fecundities are concerned, their 
mechanism cannot possibly be “basically reducible to the phenomena 
of mendelian inheritance and of a gene action’ (LERNER 1954). 
According to our results that claim is entirely unfounded, and as a 
matter of fact only rests upon a lack of adequate experiments; indeed 
reciprocal crosses, in a large gamut of external conditions, show 
sufficiently that the cytoplasm plays an important role in these 
fundamental components of fitness. But, independently from this 
statement, in our results, there is no indication for better buffering 
against temperature effects in hybrids with respect to homozygotes. 
In the hybrids the studied traits seem to be just as modifiable by 
this external factor, and their behaviour just as erratic as in pure 
strains. With respect to our basic hypothesis it means that the 
deviations of the rates of development from their nuclear Arrhenius’ 


function are just as variable and presently unpredictable as such 
deviations in pure strains. 
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At first sight, our data seem to suggest that the values of the 
Arrhenius’ constant, i.e. the temperature dependency, higher for the 
nuclear acceleration hines of the homozygotes, lower for the wild and 
hybrid strains, confirm, at least partially, the theory of genetic 
homeostasis. However this relative independence of the nuclear rate 
of the hybrids has no incidence on the degree of dependency on the 
external factors of the studied traits. This would be so if the higher 
or lower dependency on temperature of a given nucleus was suitably 
correlated and representative of the deviations of the “cytoplasmic 
rate” from the accompanying “nuclear rate” at each temperature, 
and if any fitness component of a hybrid, whether viability, develop- 
mental speed, individual size or fecundity were correlated with the 
values of the Arrhenius’ constant. This is not the case. 

Generally speaking our results lead us back to the earliest expla- 
nation of -heterosis by A. F. Smurr (1912) and G. H. Smurr (1914). 
For these authors heterosis was the result of “interaction between 
mendelfan genes and the egg-cytoplasm’’ (A. F. SHurr 1912). If one 
adds to this picture also a chromosomal interaction it seems still to 
be the most valuable approach to the problem. 

The discussion of our results with respect to the phenomenon of 
heterosis evokes directly a second problem which is touched upon by 
our results. 

In the first place all our experiments have shown that cytoplasmic 
constituents play, in first generation hybrids, a considerable role in 
the determination of characters which are all important components 
of fitness. This means in any case that the cytoplasm cannot be 
ignored on the scene of the action of natural selection. It means in 
turn that the selective value of a genotype cannot be described 
adequately without reference to a given cytoplasm, nomore than 
without reference to a given environment. It is therefore highly 
improbable that a “population” can be described in terms of “an 
integreted gene pool’ (DOBZHANSKY 1948, 1950, 1952; DOBZHANSKY 
and WALLACE 1953; WALLACE 1955; VETUKHIV 1953, 1954, 1956; 
DOBZHANSKY, PAVLOVSKY, SPASSKY and SPASSKY 1955), or, rather, 
that its desintegration by hybridization with other populations 
would be due, only, to chromosomal interactions. 

Our data have a further bearing on the problem of the hereditary 
role of the cytoplasm. The repeated backcrosses between Canton-S 
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and Swedish-B flies have demonstrated that, although the influence 
of cytoplasmic constituents is, trough reciprocal crosses, readily 
recognizable in a first generation, these reciprocal differences are 
transformed under the influence of a foreign nucleus in a given 
cytoplasm. That transformation appears in the case of the Swedish-B/ 
Canton-S backcrosses, as an “adaptation”’ of the cytoplasmic function 
to the nuclear temperature function. Indeed the replacements after 
12 generations of the Swedish-B nuclear temperature function by a 
new nuclear function toward the Canton-S nuclear temperature 
function, and vice-versa, is accompanied by the replacement of the 
Swedish cytoplasmic function by the Canton-S cytoplasmic function, 
and vice-versa. It should be noted that when we write an “adapted 
cytoplasm’’ we have no foregone conclusions about the mechanism 
or the kind of adaptation that that paraphrase conveys. 

Our experiments do not reveal the minimal number of generations 
or the minimal duration needed to transform a foreign cytoplasm by 
a given nucleus. In any case the cytoplasmic differences between 
stocks are not autonomous, i.e. they are controlled by the nucleus. 
The importance of the cytoplasm in the first generation hybrid derives 
from the fact that the action of a new nucleus upon the cytoplasm 
is not immediate. In another experiment we have shown, moreover, 
that chromosomal differences, probably limited to no more than a 
single locus, and possibly its immediate neighbourhood, is capable to 
maintain a difference in cytoplasmic behaviour. Indeed in the case 
of the +-(+)/ee(e) repeated backcrosses, it can be safely assumed 
that after 12 generations (see table 1) and with the backcrossing 
scheme used the nuclei of the +(+) and ee(e) strains, on the one 
hand, and of the +—+(e) and ee(+) on the other, are identical, 
exception made for the e gene and its immediate neighbourhood. 
For as far as our methods justify that conclusion, the (+) cytoplasm 
of the +(+) strain and the (e) cytoplasm of the ee(e) strain, 
although having been submitted to almost complete homogeneization 
of the nuclei, exception made for the +-- locus as opposed to the 
ee locus, and vice-versa, remain unchanged constituents. 

The attitude of geneticists towards the cytoplasm manifest two 
clearly discernable trends which have a tendency to move to extreme 
attitudes: or to the complete negation of a possible cytoplasmic 
influence, or to a recurrent overestimation of its importance. Judging 
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from the amount of experiments, especially on heterosis, where 
reciprocal crosses are not performed, or, if they are, are lumped 
together, most geneticists belong actually to the first category. A 
second group takes into account the true influence of the cytoplasm 
on the expression of a given organism, but goes to the opposite 
extreme and, without denying the essential hereditary importance of 
the genes, wants at any cost to charge the cytoplasm with hereditary 
properties, with autonomous hereditary constituents, for which no 
definite proof has been given. After the pioneer work of RENNER 
and KurPeRr (1921), MicHArrIs and his school have done very 
impressive work with Epilobiwm on crosses of species and geographic 
varieties with backcrosses up to 25 generations (MicHArLIs 1925, 
1935, 1939, 1943, 1950, 1951, 1953. There is a good review by 
MIicHarELIs 1954). GOLDSCHMIDT (1955) defines the problem very 
well when he wrote: “MICHAELIS reports that nearly all the properties 
of the plant are altered by cytoplasmic differences. This variability 
is much lower than that which would be produced by a series of gene 
mutations. MicHAELIs concluded that this is a definite proof that the 
cytoplasm contains autonomous genetic constituents. “GOLDSCHMIDT 
continues “T cannot concur in this conclusion. The type of effects 
studied suggests a completely unspecific cytoplasmic action. The 
cytoplasm or plasmon does not control the development of definite 
hereditary traits. It is rather a substrate for the action of the genic 
material’’. 

Whatever GOLDSCHMIDT might have meant by an ‘““unspecific’’ 
action of the cytoplasm, our data indicate that, in Drosopmla, no 
differential action of this substrate on the measured traits subsists 
without concurrent nuclear divergence, and as such they are at 
variance with MicHArEuIs’ results. Possibly there exists more widely 
divergent substrates in different species such as used by MICHAELIS. 
Maybe our Canton-S and Swedish-B strains are not even comparable 
to the geographic varieties which he studied, because of the very 
long exposure during a high number of generations to a similar 
laboratory environment. A last supposition is that the cytoplasm is 
partly made up of “permanent” hereditary components, and that it 
is further responsible for a “transitory”’ function such as that which 
we have found in our intraspecific Drosophila crosses. Whatever the 
explanation of the different behaviour of Epilobtum and Drosophila 
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cytoplasms may be, it is clear that the cytoplasmic differences be- 
tween both Drosophila strains, such as they actually are, are levelled 
out by the action of similar nuclei. They can therefore, not be called 
hereditary, not because of their apparent lack of autonomy — not a 
single process, not even the reproduction of chromosomes is ever 
autonomous in living being — but because the processes involved in the 
reproduction of these cytoplasmic differences is to a much lesser 
degree constant in the course of time than those of chromosomal 
differences. 

The situation in Drosophila can be-adequately summarized by 
MORGAN'’s sentence: “The cytoplasm may be ignored genetically” if it 
is well understood that it cannot be ignored at all in genetical ex- 
periments. In such experiments, and more explicitely in the first gener- 
ation, the cytoplasm plays the role of a “primary environment” 
which can lead to quantitatively different phenotypic expressions. 
There is no reason why it should not be able to switch genetic 
potentialities in qualitatively different pathways, exactly like the 
“secondary environment’ can switch the dominance relationship 
between alleles. Such a role the cytoplasm seems to play in the 
realization of mating types in Paramecium (SONNEBORN 1954; see 
also NANNEY 1957). We can, therefore, not agree with SONNEBORN 
when he insists, on the basis of his results, on the inclusion of the 
cytoplasm “as genetic material and as a part of the genetic system” 
(SONNEBORN 1954). 

Although our results disclose new ways of tackling the problem of 
heterosis, of cytoplasmic action, and of quantitative inheritance in 
general, a number of important questions are left open after our 
investigation, or received only a partial answer. 

À first question concerns the relationship between nucleo-cyto- 
plasmic rates and adaptation. SPAAS and Heuts (1958) find in 
Salmonidae precise temperature conditions where given forms ex- 
hibit a viability of 100% or nearly so. These temperatures to which 
the eggs are clearly best adapted correspond to points were develop- 
mental rates cross the hypothetical nuclear acceleration line. In our 
experiments the disposition of developmental rates around the 
nuclear rates are much irregular, and very frequently do not cross 
the nuclear rates at any point within the temperature range studied. 
It is clear that, on the one hand, the number of experimental temper- 
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atures in our experiments was insufficient to get an idea of the real 
shape of the cytoplasmic acceleration curve, and, on the other hand, 
that, in contrast to SPAAs and Heuts’ experimental conditions 
competition plays an important role. This competition certainly has 
an influence on the developmental rate, and may be such as to shift 
the developmental rate curve, without, or without as much dis- 
placement of the nuclear curve i.e. without a change of the parameters 
of the Arrhenius’ function. It is further to be expected that an 
identical population density at different temperatures does not result 
in an identical degree of competition at each of them. The degree of 
competition will depend in any case on the different metabolic rates 
under these conditions. Repetition of our experiments at more 
temperatures and without any competition is indicated. 

A second point is that as yet our data are insufficient to explain 
the greater steepness of the nuclear acceleration line in relation to the 
increasing homozygosity of the genome, nor to explain the relatively 
low value of the Arrhenius’ constant of the Inbred strain in relation 
to the Canton-S strain from which it was derived. It is likely that the 
answer to that question lies in the relation between the increase of the 
Arrhenius’ constant values and the increase in steepness of the 
cytoplasmic curve in function of the number of inbreeding generations. 
Maybe a solution can be given to that problem in the study at the 
points of view of development speed and viability, generation after 
generation, and at diverse temperatures, of diverse laboratory and 
wild strains submitted to inbreeding. Lastly let us point out that if, 
in the hybrid, the fact of interaction between the genomes of the 
parental strains is well established, nothing is yet known about the 
mechanism of that interaction, although some indications for a 
“dominance”’ of one genome on another in the establishment of the 
nuclear rate of the hybrid exist. 

A last important question concerns the causal relationship between 
the different traits studied (fig. 14). Our experiments prove the 
existence of several correlations between these traits, but to describe 
these connections in terms of cause and effect is a rather delicate task 
for which, at this phase of our work, a number of new hypotheses 
have to be advanced. We have supposed that the rate of development 
is to be considered as a resultant of two independent and specifically 
different thermosensitive systems whose location we have traced back 
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to the nucleus and to the cytoplasm. Viability, according to our 
assumption, is directly related to the deviations of both hypothetical 
rates from each other at each temperature. As such viability is 
related to the developmental rate but this relationship is spurious, 
both traits being independently caused by the deviations mentionned, 
for which we have used the developmental rate as a technical 
measurestick. Egg viability itself modifies the environment by 
causing variable population densities, and, in a competitive en- 
vironment, a variable degree of competition, which in turn influences 
the larval survival. The final viability is expected to influence in a 
negative way the average size of the survivors; to such an effect we 
trace back a negative correlation between size and viability. It is 
conceivable that this secondary environmental influence conceals some 
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Fig. 14. Causal links between the different traits studied. Explanation in text. 


direct relationship between final size — independent from the develop- 
mental speed, as defined for viability — and the nucleo-cytoplasmic 
rate deviations. New experiments without competition might give an 
answer to this open question. A similar supposition can be made 
concerning fecundity, were variation due to zygotic causes might 
have been concealed by the effects of food conditions. Such an 
assumption is supported by the fact of a positive correlation between 
fecundity and viability, unless it is supposed that this correlation is 
due to some positive effect of population density on fecundity. If 
fecundity is directly related to the nucleo-cytoplasmic rate relation- 
ship a positive relationship with developmental rate as well as with 
nucleo-cytoplasmic rate deviations is expected. We consider, however, 
our fecundity measurements as too incomplete to push, for the 
moment at least, our argumentation any further in this direction. 
A last remark concerns the relation between fecundity and size. 
This last trait appears to be largely dependent on external conditions, 
chiefly those provoked by a higher or lower viability, resulting in a 
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negative correlation between both values. The residual variation of 
size around the regression line of size on viability proved however 
not to be without significance, in as much that this variation was 
found to be negatively and significantly correlated with a similar 
deviation of the fecundities from the size on fecundity regression 
line, while there is no relationship between the first deviations with 
the fecundity deviations from a fecundity on size regression line. In 
other words a residual size variation which is not concealed by the 
influence of competition proves to be related inversely to similar 
variations in fecundity around a regression line, assuming a dependen- 
cy of size on fecundity such as to correlate a plus variant in fecundity 
to a minus variant for size. This rather unexpected result would not 
seem as unorthodox in the light of a further hypothesis, namely that 
the primary genetic trait of a given zygote would be the potentiality 
to devide a given number of times, with a given speed. This potenti- 
ality could be unequally and complementarily attributed to the germ 
track and to the somatic cells. This hypothesis would require, how- 
ever, that a greater size, for as much as it is related to lower fecundity 
would be the consequence of a higher cell number, which remains to 
be demonstrated. 


SUMMARY 


L. Viability, development time and development rate, size and fe- 
cundity of diverse Drosophila melanogaster homozygotes and hybrids 
strains, have been studied in a large gamut of experimental conditions, 
including particularly variations in temperature and systematic study 
of reciprocal crosses. 

2. The relativity of the traits of homozygotes and hybrids is demon- 
strated in relation 1) to the ages and development times of the pa- 
rental strains, 2) to the variation of the temperature conditions, 3) 
to the direction of reciprocal crosses. 

3. It is shown that at each temperature the mortality of a particular 
nucleo-cytoplasmic combination, of which the cytoplasm and the 
nucleus are of the same origin, is determined by the interaction of 
two different thermo-reactive rate systems, the first being a property 
of the nucleus, the second a property of the cytoplasm. These thermo- 
reactive developmental rates, are graphically represented and theo- 
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retically defined as “nuclear acceleration line’ and “cytoplasmic 
acceleration curve”. 

4. The mortality of a hybrid is shown to be determined by the 
same relationship, the nuclear rate being a function of the interaction 
between the parental genomes and the cytoplasmic rate a function of 
the development speed of the maternal parent. 

5. Nuclei which are genetically different, at least for a single 
mutation possess a different cytoplasm. This cytoplasmic differences 
are shown to be dependent, however on the nucleus, and therefore 
the cytoplasm cannot be considered as an hereditary constituent. 

6. Size and fecundity are, in opposite ways, dependent on the 
viability, or more precisely on the environmental conditions created 
by a higher or lower viability. They are also, again in opposite ways, 
dependent on zygotic factors, i.e. the relationship between nuclear 
and cytoplasmic rate. In the experimental conditions used size is 
more dependent than fecundity on viability and on the zygotic factors. 

7. The nuclear rate can be adequately expressed by the Arrhenius’ 
constant. This constant is different for males and females of the same 
strain. For a particular strain the value of the Arrhenius’ constant 
does not allow any inferences concerning the quantitative and 
homeostatic properties of that strain. 

8. The results are discussed in relation to the diverse theories of 
heterosis, and with respect to the relative importance of the cytoplasm 
versus the nucleus in heredity. Finally an hypothetical scheme of 
causal relationships between the different quantitative traits studied 
is presented and discussed. Some repercussions of the results on future 
research in quantitative genetics are considered. 
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HISTORICAL NOTES, III 


M. J. Sirks (Genetisch Instituut R.U., Haren (Gron.) Netherlands) 
HAECKEL IN 1866 ON THE ROLES OF PLASMS AND NUCLEUS 


(HAECKEL, ERNST, 1866. Generelle Morphologie der Organismen 
(Berlin, Reimer, 1866), 1 : 287-288). 


„Hier scheinen nun die beiden Funktionen der Erblichkeit und der Anpassung 
bei den kernlosen Cytoden noch nicht auf differente Substanzen vertheilt zusein, 
sondern der gesamten homogenen Materie des Plasma zu inhärieren, während 
dieselben bei den kernführenden Zellen in der Weise auf die beiden heterogenen 
activen Substanzen vertheilt sind, dass der innere Kern die Vererbung der 
erblichen Charaktere, das äussere Plasma dagegen die Anpassung, die Accomo- 
dation oder Adaptation an die Verhältnisse der Aussenwelt zu besorgen hat … 
Wenn wir demgemäss das Plasma vorzugsweise als den nutritiven, den Nucleus 
dagegen vorzugsweise als den reproduktiven Bestandtheil der Zelle ansehen 
können, so werden wir mit Recht den Kern der Zellen als das hauptsächliche 
Organ der Vererbung, das Plasma als das hauptsächliche Organ der Anpassung 
betrachten können”’. 

Translation into english: 

‚In the anucleate Cytodes both functions, that of inheritance and that of 
adaptation do not seem to be distributed among differentiated substances, while 
they both are performed by the homogeneous substance of the plasm alone; in 
nucleate cells however these functions are separately entrusted to both hetero- 
geneous active substances in this way, that the inner nucleus is charged with the 
inheritance of hereditary characters, the exterior plasm takes care of the acco- 
modation or adaptation to the conditions of environment … So, if we consider the 
plasm chiefly the nutritive part of the cell, the nucleus however the reproductive 
part, then the nucleus of the cells may be justly regarded as the principal organ 
of inheritance, the plasm as the principal organ of adaptation”’. 
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